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ABSTRACT 

We report on the local environment, radial velocity, and star formation history of [O ll] -selected emission- 
line objects in the z ~ 0.4 galaxy cluster Abell 851, using optical spectra obtained with Keck I/LRIS and 
Keck II/DEIMOS. A large fraction 55%) of cluster [Oil] emitters show strong Balmer absorption lines 
(> 4 A in equivalent width at i\5). These e(a)-type spectra have been attributed to dusty starburst galaxies 
by Poggianti & Wu (2000), an interpretation supported by our Balmer-decrement-derived reddening measure- 
ments, which show a high frequency of very reddened [E{B-V) > 0.5] [O ll] emitters. Our spectral modeling 
requires starburst ages < 1 Gyr, which is shorter than the cluster crossing timescale. We argue that this star- 
burst phase occurs during cluster infall based on the radial velocity distribution of the cluster [O ll] emitters, 
which present a deficit of systems near the cluster systemic velocity compared to a virialized population (or a 
backsplash population). The spatial segregation of some redshifted and blueshifted groups strongly indicates 
that the accretion was recent. Throughout the cluster, the fraction in [O ll] -emitting galaxies is a strong function 
of the local galaxy density, being suppressed in dense environments. Our analysis supports previous sugges- 
tions that dusty starburst galaxies arise at the expense of continuously star-forming gas-rich spiral galaxies [i.e., 
the e(c)-type; Dressier et al. (1999)]. In addition, we describe a fainter [Oil] -emitting population, comprised 
largely of dwarf galaxies (Martin et al. 2000) and find an even stronger suppression of [O ll] emission in high 
density environments among this subsample, indicative of more effective destruction by harassment and/or gas- 
stripping. Comparison to previous morphological studies, limited to the central region of Abell 851, suggests 
that galaxy-galaxy interactions may trigger the starbursts. The high e(a) galaxy fraction in Abell 85 1 compared 
to that in the field, however, suggests that some cluster-specific mechanism, likely related to the dynamical 
assembly of the cluster, also contributes to the high number of starbursts. 

Subject headings: galaxy cluster: general — galaxies: clusters: individual(Abell 85 1 ;CL0939-i-47 13) — galax- 
ies: dwarf — galaxies: evolution — galaxies: starburst — cosmology: observations 



1. INTRODUCTION 

Local environment has a strong effect on the star-forming 
activity and morphology of galaxies — the observation first 
recognized in the cores of rich clusters (Dressier 1980; 
Butcher & Oemler 1984). There is mounting evidence in sup- 
port for the Butcher-Oemler effect (Butcher & Oemler 1984) 
as a cosmic process, not specific to galaxy clusters, in which 
the field-to-cluster transformation of galaxies occurs upon the 
accretion of field galaxies to cluster central regions (Dressier 
et al. 1994; Couch et al. 1994; Abraham et al. 1996; ElUng- 
son et al. 2001; Tran et al. 2004). In fact, recent studies 
of cluster galaxies have shown that mechanisms responsible 
for their transformation must be in effect well beyond cluster 
virial radii, well before the galaxies at infall reach the cen- 
tral regions of clusters (e.g., Balogh et al. 1997; Treu et al. 
2003; Rines et al. 2005; Moran et al. 2005; Pimbblet et al. 
2006). Since cluster environment samples a wide range in lo- 
cal galaxy density and relative velocity of cluster members, 
a dominant physical process driving galaxy evolution may 
vary over an infall period. Thus most recent efforts in clus- 
ter galaxy surveys have focused on studying galaxies in infall 
regions where the preprocessing of cluster galaxies is already 
in effect. 

Although the studies spanning a few decades have estab- 
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lished such empirical relations as the morphology-density re- 
lation (Dressier 1980) and the Butcher-Oemler effect (Butcher 
& Oemler 1984) for giant luminous cluster galaxies, the un- 
derlying physical processes that are responsible for the trans- 
formation of their properties and therefore the aforementioned 
empirical relations have yet to be well constrained to date. 
The sign of the differential evolution between dwarf and gi- 
ant galaxies is even more obscure, since few wide-field sur- 
veys have sampled faint galaxies at high redshift. Dwarf 
and large spiral galaxies are expected to respond differently 
to cluster tidal fields due to their different mass concentra- 
tions (Moore et al. 1999). In predicting a steep halo mass 
function, the currently standard hierarchical galaxy formation 
models (e.g., Kauffmann et al. 1993; Cole et al. 1994) pro- 
duce more dwarf galaxies than actually observed (e.g., Klypin 
et al. 1999; Somerville 2002). Galaxy clusters host dense en- 
vironments that are hostile to fragile dwarf galaxies; therefore 
the apparent faint-end steepening of cluster galaxy luminosity 
function presents a challenge to dwarf formation models (see, 
e.g., Popesso et al. 2005, and references therein). 

Redshift surveys of local universe and A^-body simulations 
show that large-scale filaments are a natural product of grav- 
itationally driven structure formation (e.g.. Bond et al. 1996; 
Rines et al. 2001; Pimbblet et al. 2004). Galaxy clusters are 
expected to form at the intersections of filamentary structures, 
yet there has been few studies on the direct detection of fila- 
ments extending out from isolated clusters at high redshifts 
(e.g., Ebeling et al. 2004). A rich, moderately distant (z ~ 
0.4069) galaxy cluster Abell 851 [CL0939H-4713; 9''42™56?2, 
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46°59'12"0 (J2000)] offers evidence of filaments, observed 
via the photometric redshift technique (Kodama et al. 2001, 
hereafter KOI). Since few clusters at this epoch have been 
studied well beyond their virial radii, it remains unclear 
whether the observed structures are actually filaments of cos- 
mological origin or are related to on-going merger event (De 
Filippis et al. 2003). In either case, Abell 851 presents a 
unique opportunity to study cluster galaxy evolution in a dy- 
namically active environment on the large scales. 

Before the identification of the filaments, Abell 85 1 had al- 
ready drawn attention due to the extreme richness (e.g., Gunn 
et al. 1986) and high abundance of blue, star-forming galax- 
ies; the cluster contains a larger fraction in young emission- 
line galaxies compared to clusters at similar redshift (Lotz et 
al. 2003). A spectroscopic survey of r < 22.5 galaxies re- 
vealed a high fraction in post-starburst galaxies, whose spec- 
tra are dominated by Balmer absorption lines of intermediate- 
age stars (Dressier et al. 1999, hereafter D99). The timescales 
associated with such galaxies impUed the suppression of star- 
forming activity within the past ~ 1.5 Gyr. Since the clus- 
ter crossing time for infalling galaxies is at least a few Gyr 
(Treu et al. 2003) and definitely longer than the post-starburst 
timescales, the high abundance of post-starburst galaxies sug- 
gested the significance of cluster-induced starburst in the 
transformation of cluster galaxies. Although the scarcity of 
star-forming activity in very dense environments within the 
central regions of clusters has been universally observed and 
confirmed, there has been a debate as to the existence of en- 
hanced starburst phase before the suppression of star-forming 
activity of infaUing cluster galaxies kicks in (Dressier et al. 
1999; Balogh et al. 1999; Tran et al. 2003, 2004). In post- 
starburst galaxies, we see remnants of starbursts, but detecting 
the signs of starbursts is generally complicated by the need for 
constraining the star formation history (SFH) of galaxies — 
the process requiring more than a simple diagnostics of star- 
forming activity. 

Martin et al. (2000, hereafter MLFOO) carried out narrow- 
band and broadband photometry of Abell 851 using May- 
all 4-meter telescope at the Kitt Peak National Observatory, 
to detect faint emission-line galaxies using their redshifted 
[Oii]A3727 emission Une. The narrowband search yielded 
371 emission-line cluster member candidates. The complete- 
ness limit of the MLFOO survey was g ~ 24.5, about ten times 
deeper than the D99 spectroscopic survey. The locus in the 
{g — i)-i color-magnitude diagram and their small isophotal 
area suggested their dwarf identity for a significant fraction of 
the [O II] emission candidates. In a spectroscopic follow-up 
program of the MLFOO survey, we have obtained an exten- 
sive set of spectra of the |0 ll|-emitting candidates using the 
Keck I Low Resolution Imaging Spectrometer (LRIS; Oke et 
al. 1995) and Keck II Deep Imaging Multi-Object Spectro- 
graph (DEIMOS; Faber et al. 2003). 

As the second in a series of papers drawing on the Keck 
spectroscopic sample of the [Oil] emitters, we present the 
analysis of environment, kinematics properties, and SFH of 
infaUing faint star-forming galaxies in a highly dynamically 
active galaxy cluster at z ~ 0.4. The observations, data, and 
the spectroscopic confirmation of MLFOO [O ll] emitters have 
been reported in Sato & Martin (2006, hereafter Paper I), to 
which readers are referred for the detailed information on the 
basic characteristics of our spectroscopic sample. This paper 
starts with the sections entirely dedicated to the effect of local 
envirormient(§ 2) and the kinematic properties (§ 3) of cluster 
[O II] emitters. In § 4 we present the cluster [O ll] emitters in 



terms of their spectral types derived from the [O ll] emission 
and H(5 absorption equivalent widths, which will be analyzed 
in view of of photoionization and stellar population synthesis 
models in § 5. We then discuss (§ 6) and sunmiarize (§ 7) our 
results. 

Throughout this paper, we use the standard ACDM cosmol- 
ogy, (f]„,riA) = (0.3,0.7), with Ho = 70 km s-i Mpc"'. At the 
distance of Abell 851 (z = 0.4069), 1" on the sky corresponds 
to a projected physical distance of 5 .43 kpc, and the lookback 
time to the galaxy cluster is 4.3 Gyr. 

2. EFFECTS OF LOCAL ENVIRONMENT 

The star-forming property of a galaxy has been shown to 
depend strongly on its local environment, often quantified by 
the projected number density of galaxies near the apparent 
position of the galaxy in the sky. For Abell 851 in particu- 
lar, KOI have reported a so-called threshold effect — galaxy 
V-I colors abruptly become redder above a threshold den- 
sity, ~ 50 galaxies Mpc"^ (in the cosmology assumed in this 
paper), which roughly corresponds to the typical local galaxy 
density within the substructures in surrounding filaments.^ 
Robust determinations of cluster memberships are essential 
for quantifying local environment of cluster galaxies. We ob- 
tained the photometry and photometric redshifts of the objects 
studied in KOI, kindly provided by T. Kodama. The KOI sur- 
vey limiting magnitude was roughly comparable to our sur- 
vey: /(AB) ~ 23.4 (KOI) and / ~ 23 (MLFOO); this is quite 
deep, as KOI estimates My to be about /= 19.4 in this cluster. 

2.1. Local Projected Density of Galaxies 

The local environment of a galaxy is quantified by a pro- 
jected surface galaxy number density, S, using ten nearest 
cluster members defined by KOI via the photometric redshift 
method. The value is computed for each galaxy in our sur- 
vey as follows: (1) Using the RA and DEC of a galaxy, we 
find ten KOI confirmed cluster members nearest to that posi- 
tion; (2) compute the distance between the farthest and closest 
members to that reference coordinates in the ten objects just 
found; and (3) take that as the radius of a circular projected 
area by which the number of galaxies within the area (i.e., 
always 10) is divided. The only difference from the conven- 
tional method is that our ten-nearest neighbor density may or 
may not include the galaxy for which the value is computed. 
This is because explicit spatial cross-correlation was not made 
between the MLFOO and KOI objects, but the deviation does 
not significantly alter the conclusions drawn from our study. 
KOI estimated the field contamination to the projected density 
of cluster galaxies in their survey to be 2.29 ± 0.06 arcmin"^. 
We statistically subtracted this from our computation of local 
galaxy density described above. In the nearest-neighbor ap- 
proach, the effect of Poissonian uncertainty in the expected 
number of field galaxies becomes dominant for smaller S, 
since the area enclosing ten cluster members contains many 
more non-member objects. A Icr Poisson uncertainty in the 
field contamination becomes comparable to the local density 
estimate for a cluster member at logS ~ 0.3. The system- 
atic removal of field contamination yields unphysical value of 
S < Mpc~^ for some galaxies. Where required, we gener- 
ally assigned logS = -2 as the upper limit for such galaxies. 

' We caution that local galaxy densities computed in such a method used 
in our study may not be directly compared to those in other surveys; a pro- 
jected density estimate derived from a nearest-neighbor approach is highly 
sensitive to the specific ways in which foreground and background galaxies 
are removed from a cluster sample. 
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Fig. 1 . — (Top) The [O II] emitter fractions (filled circles) as a function of 
local galaxy density S. The fractions of [O II] emitters in g — i < 1 galaxies 
(open squares) are also shown. The fractions are computed from the number 
of spectroscopically-confirmed [O ii] emitters (multiplied by appropriate cor- 
rection factors given in Paper I to account for the selection bias) normalized 
to the nvmiber of KOI cluster members at each density bin. Dotted horizon- 
tal lines indicate the best estimates and la uncertainties for the fractions of 
galaxies with Wo([0 II]) < -10 A in cluster and field samples at z ~ 0.4 com- 
puted by Nakata et al. (2005). (Bottom) Fractional abundance as a function of 
local galaxy density for the spectroscopically-confirmed cluster [O II] emit- 
ters corrected for the target selection bias (N ~ 175; solid histogram) and 
the KOI cluster members (N = 808; dashed histogram). The fractional abun- 
dance is normalized to the total number of objects (N) within each sample to 
facilitate the comparison of density distributions. The vertical gray dashed 
line indicates the threshold density, logS ~ 1.7, over which V—I colors of 
galaxies become redder as reported by KOI. 

2.2. [Ou] Emitter Fraction 

We start with a brief summary of the basic properties of the 
cluster [O ll] -emitting sample presented in Paper I. Our clus- 
ter [0 11] emitters for this spectroscopic follow-up were drawn 
from 371 narrowband-selected objects with their candidate 
[O II] emission detected at S/N > 3 level in MLFOO. The nar- 
rowband survey was roughly complete to the observed equiv- 
alent width of < -11 A at m5i29(AB) ~ 24.^ The preference 
was given to the objects with bluer ; colors (i.e., g-i< 1), 
so the sample is more complete for blue objects. Among 
the 212 objects selected for spectroscopic follow-up, 111 ac- 
tually showed [Oil] emission and had redshifts within the 
range z = 0.395-0.420, constituting our cluster [O ll]-emitting 
sample. This final sample spans a wide range in luminosity, 
18.5 < i < 25.5, where ; ~ 23 corresponds to the complete- 
ness limit of the MLFOO narrowband survey. 

Fig. 1 shows the distributions of cluster [O ll] emitters and 
the KOI cluster members within the survey area of MLFOO. 
Note that the histogram for the spectroscopically-confirmed 
cluster [On] emitters are corrected for the target selection 
bias by \he g — i color, i.e., we multiply the number actually 
observed by the the correction factor for each g — i color sub- 
sample; see Paper I for the details. A Kolmogorov-Smirnov 
(K-S) test of the distributions of local galaxy density for the 
spectroscopically-confirmed cluster [O ll] emitters and more 
general, photo-z-selected cluster population from KOI yields 
a very small probability (-C 1%) that they share a common 
parent distribution. This indicates the [O ll] -emitting popula- 
tion resides in different environments than the general clus- 

* An equivalent width of an emission line is designated a negative sign; 
see Paper I for the sign and other notational conventions used in this paper. 
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Fig. 2. — The [On] emitter fraction of KOI photo-z cluster members as a 
function of projected radius from the cluster center The center of the clus- 
ter was taken to be at 9'"43°"03?3, 46°59'26;'5 (J2000), corresponding to 
the brightest point in X-ray surface brightness from De Filippis et al. (2003). 
Dotted horizontal hnes indicate the best estimates and la uncertainties for 
the fractions of galaxies with Wo([0 II]) < -10 A in cluster and field samples 
at z ~ 0.4 computed by Nakata et al. (2005). At each radial bin, the frac- 
tion is corrected for MLFOO emission-line candidates not spectroscopicaUy 
observed, in addition to the usual target selection correction by g-i colors 
(Paper I). 

ter galaxy population. A relatively higher fraction of clus- 
ter [0 11] emitters is found at low densities. The emission- 
line galaxy fraction has been shown observationally to de- 
pend strongly on local environments, where a smaller fraction 
is observed in dense environments (e.g., Balogh et al. 2004b). 
Furthermore, such an emission-line galaxy fraction appears to 
evolve with redshift (Kodama et al. 2004; Nakata et al. 2005). 

At the lowest density bin, the fraction in [On] emitters 
matches well the field [Oil] emitter fraction at z ^ 0.4 of 
~ 70% estimated by Nakata et al. (2005) from their sam- 
ple drawn from the Canadian Network for Observational 
Cosmology 1 (CNOCl) cluster survey; the [On] equivalent 
width cut for their [O ii]-emitting population is similar to ours 
[Wo([0 II]) < -10 A]. As discussed in Paper I, MLFOO could 
have detected only ~ 75% of virialized cluster |0 II] emitters 
with their narrowband filter, assuming a Gaussian velocity 
distribution for viriaUzed cluster members. In Fig. 1, we did 
not statistically correct for this population which could have 
been missed by the MLFOO narrowband survey. Although we 
do not expect the MLFOO survey area (4.4 Mpc x 4.4 Mpc) to 
be large enough for sampling a true field population and ap- 
pUed a crude correction for target selection bias, the sampling 
at low-density region appears to provide us with a reasonable 
[O II] -emitting sample in field-like environments, since their 
observed fraction is consistent with the independent measure- 
ment of field [O ll] emitter fraction. 

The [O ll] emitter fraction in the densest environment in our 
sample, ~ 20%, is also roughly consistent with the observed 
[On] emitter fraction in z = 0.4 clusters (Nakata et al. 2005). 
In relation to the KO 1 threshold density at which the V-I col- 
ors of cluster galaxies have been shown to become red quite 
abruptly, the presence of [O ii] emitters appears already sup- 
pressed at slightly lower densities, and their fraction remains 
roughly constant at the lowest level toward denser environ- 
ments. Fig. 2 shows the [O ll] emitter fraction as a function 
of projected radius from the cluster center. The radial trend 
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Fig. 3. — Cumulative local density distributions of the cluster [O II] emit- 
ters with various D99 spectral types: e(b) (open squares), e(c) (open circles), 
and e(a) (open triangles). The distribution is normaUzed to zero at the field 
value (i.e., at the lowest density bin) for each spectral type. The samples are 
binned with constant logS = 0.5 intervals, where the symbols in the plot in- 
dicates their central values. The spectral types are assigned using the best 
estimates as well as the uncertainties in [O II] and US measurements (see 
§ 4.2 for detail). The fractions are corrected for the target selection bias by 
appropriate correction factors (Paper I). In the electronic edition, the relative 
abundances of the [O II] cluster emitters in different D99 spectral types are 
also indicated by the filled curves: from the top, they are e(b) (blue), e(c) 
(green), e(a) (red), k (yellow), and k-l-a/a-Hk (magenta). The colors are half- 
toned for better contrast. [See the electronic edition of the paper for a color 
version of this figure.} 

is expected in view of Fig. 1, since the local galaxy density is 
generally anti-correlated with the cluster radius. 

The distribution of ^ - / < 1 galaxies among the [O ll] emit- 
ters shows a modest peak near the KOI threshold density. It 
is hard to see whether this peaking is directly related to the 
threshold effect, but it could indicate the presence of star- 
bursts as galaxies enter local environments typical of cluster 
substructures. In Fig. 3, we see the cumulative local density 
distributions of the cluster [Oil] emitters with various D99 
spectral types. According to tests, the cumulative distribu- 
tions among e(b), e(c), and e(a) types are different at very high 
confidence levels 99%). The details of the D99 classifica- 
tion scheme will be discussed in § 4.2, but for the current dis- 
cussion, we note that e(b), e(a), and k-i-a/a-i-k galaxies are all 
associated with starbursts, whereas e(c) and k galaxies are not. 
Fig. 3 may provide additional evidence for a higher incidence 
of starbursts in the local environments near the KOI threshold 
of log E ~ 1.7. We also note a striking drop in the [O ll] emit- 
ter fraction at log S « 1.3, not too far from the KOI threshold 
density. Although a very strong suppression of star-forming 
activity near the KOI threshold could be real, the [O ll] emitter 
and g-i < 1 galaxy fraction are already well below the field 
values at logS « 0.5; therefore, some mechanism of quench- 
ing star formation could be in effect below the KOI threshold 
density as well. 

2.3. [Oil] & Ha Equivalent Width Distribution 

Although the star-forming galaxy fraction strongly depends 
on local environments, in previous studies the amount or 
strength of star-forming activity within each galaxy, probed 
via such emission-line diagnostics as Ha equivalent width, 
does not appear to depend strongly on local environments 
(Balogh et al. 2004a,b; Kodama et al. 2004; Rines et al. 2005). 
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Fig. 4. — Rest-frame [O II]A3727 (bottom) and Hq (top) equivalent widths 
of cluster [O II] emitters as a function of their local galaxy density. Open and 
filled circles are H II region-like galaxies and AGN-like objects, respectively; 
crosses are the objects whose classification by the line ratio diagnostics could 
not be carried out. See Paper I for the method used for these identifications. 
For the data points with log S = -2, small negative random numbers are added 
to their densities to avoid excessive overlapping. 

In Fig. 4 we see how our cluster [O ll]-emitting sample com- 
pares. For the wide dynamic range in local galaxy densi- 
ties, almost five decades in S, the strength of [O ll] and Ha 
equivalent widths shows only a mild decrease at high densi- 
ties. The lack of dependence on local galaxy densities be- 
comes more remarkable in view of Fig. 5, where luminous 
[O II] emitters — many of which are found to be AGN-like (Pa- 
per I) — have weaker [O ll] and Ha equivalent widths. The 
high abundance of luminous AGN-like [O ll] emitters at high 
densities significantly accounts for the apparent trend. 

Relative independence of the strength of emission equiv- 
alent widths of star-forming galaxies on their local galaxy 
density, combined with the strong trend of decreasing [O ll] 
emitter fraction with local galaxy density, has given rise to an 
interpretation that the mechanisms that halt star-forming ac- 
tivity in these galaxies either work on very short timescales 
and/or are sufficiently rare that we are not seeing an abun- 
dance of galaxies in transition from actively star-forming to 
passively-evolving states. Starbursts followed by quench- 
ing of star-formation offers an attractive solution in this re- 
gard. A fraction of total galaxies in the starburst or post- 
starburst phase could offer an insight, yet the existing surveys 
have shown considerable disagreements as to the frequency of 
starburst and/or post-starburst galaxies (Dressier et al. 1999; 
Balogh et al. 1999; Tran et al. 2003, 2004). Even with a 
robust estimate of post-starburst fraction, the role of post- 
starburst phase in galaxy evolution is still obscure. Nonethe- 
less our cluster [O ll] -emitting sample offers no contradicting 
evidence to either of the above interpretations. 

2.4. Luminosity Distribution 

In Fig. 5 we see the distribution of cluster [O ll] emitters 
in terms of their luminosity and local environment. Interest- 
ingly, it shows that among our cluster [O ll] emitters AGN- 
like galaxies are luminous and mostly reside in dense envi- 
ronments. However, we must caution that within the inter- 
mediate luminosity range there are a number of galaxies for 
which we could not carry out line-ratio diagnostics for AGN 
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Fig. 5. — The (■ magnitudes of cluster [On] emitters as a function of their 
local galaxy density. Symbols are as in Fig. 4. Grids of the ;-S plane defines 
the subsamples of cluster [On] emitters to be coadded (§ 5.2). The magni- 
tude cuts are made at = 21 and 23, and the density cuts are at log S = and 
1.6. Objects are divided into nine subsamples by their log S: low-density (L), 
transitional (T), and high-density (H) environments; and by their luminosity: 
bright (B), medium (M), and faint (F) in term of i magnitude. 

activity (Paper I), indicated by crosses in the figure. Also, 
none of thie cluster [O llj emitters are spatially coincident with 
the XMM-Newton X-ray point sources detected at the level 
of Lx(bol) > 1.6 X 10*3 erg s"' by De FiUppis et al. (2003). 
The X-ray observation does not exclude the possibility of 
these galaxies hosting nuclear activities, however. In Fig. 6 
we see the AGN-like galaxies in our sample have g-i > 1, 
putting them in the reddest end of the distribution of SDSS red 
AGNs, which generally appear X-ray weak as well (Brandt et 
al. 2004). If the observed trend of AGN-like galaxy distribu- 
tion is real, it would point to the significant role of optically- 
luminous AGN-like galaxies in dense environments, whose 
[O II] emission equivalent widths are generally smaller than 
those of more general star-forming galaxy population (Pa- 
per I). As mentioned in Paper I, if AGN feedback acts to sup- 
press star formation in luminous (or massive) galaxies, that 
offers a possible explanation for the downsizing effect, which 
refers to observational evidence for the redshift evolution of 
characteristic galaxy luminosity above which star-forming ac- 
tivities are suppressed (Cowie et al. 1996). 

It is also notable that we see relatively few low-luminosity 
galaxies in dense environments, in the region marked as HF 
in Fig. 5. Although the range i > 23 is below the complete- 
ness limit of narrowband survey, the depletion of faint galax- 
ies at high densities is remarkable because of the abundance 
of galaxies with similar luminosities in lower-density environ- 
ments. The relative depletion of faint galaxies in high-density 
environments here simply means the lack of star-forming pop- 
ulation, and can result from suppressed star-forming activity 
in faint galaxies, which simply fade out of our survey limit, or 
their true absence, possibly due to the destruction of low sur- 
face brightness dwarf galaxies within a strong cluster grav- 
itational potential or by a series of impulsive gravitational 
encounters with other galaxies — so-called galaxy harassment 
(Moore et al. 1999). The j'-band magnitude of our [O il] emit- 
ters samples their rest-frame V-band light, which itself de- 
pends on the amount of on-going star-forming activity. Since 
our survey is insensitive to the non-star-forming population, 



we cannot rigorously exclude or favor either of the two possi- 
bilities. 

In Paper I we showed the strength of [O ii] equivalent width 
shows only a mild decreasing trend with luminosity. In dense 
environments, the absence of strong [O ii] emitters is partially 
a consequence of the lack of faint [O ll] emitters and the pres- 
ence of optically-luminous AGN-like galaxies with moderate 
to weak [On] equivalent widths. This environmental trend 
of [O ll]-emitting population makes stronger the claim that 
the amount of star-forming activity in star-forming galaxies 
does not vary much in different local environments, while the 
star-forming galaxy fraction is a strong function of local envi- 
ronment. The luminosity segregation of galaxy populations in 
dense environments further suggests the importance of differ- 
ential evolution of galaxies based on their intrinsic properties. 

3. KINEMATICS OF THE CLUSTER [O n] EMITTERS 

Studying the connections between galaxies and their host 
environments in galaxy clusters helps to understand the trans- 
formation of galaxies in cosmological context. Cluster galaxy 
populations in relation to projected cluster-centric radii and 
galaxy local number densities are relatively easier to quantify 
and therefore have been well studied (e.g., Oemler 1974; Mel- 
nick & Sargent 1977; Treu et al. 2003). The relations between 
cluster galaxy populations and their kinematics, however, are 
nontrivial since cluster galaxies are not necessarily in dynam- 
ical equilibrium on the scale of an entire cluster. Whereas 
projected cluster-centric radii set lower limits on the true dis- 
tances of cluster members from the cluster dynamical center, 
observed line-of-sight velocities say little about their relative 
locations to the cluster, not to mention obtaining a number 
of redshifts requires observationally expensive spectroscopic 
surveys. As the connections between the morphology-radius 
and morphology-density relations and the hierarchical picture 
of structure formation become clearer, it is of great interest to 
better constrain the relations between the dynamical states of 
clusters and their effect on constituent cluster galaxies. 

In this section, we analyze the kinematic states of 
cluster [On] emitters. This population is compared to 
the spectroscopicaUy-confirmed early-type galaxy population 
from the D99 survey to shed more light on the connections 
between the dynamical states and the star-forming property 
of cluster [O llj emitters. Although the size of D99 spectro- 
scopic sample is small, = 31 for k-type galaxies (see § 4.1 
for definition of spectral types), comparing the [O ii] emitters 
to the galaxies composed of much older stellar population ob- 
served in the same cluster illuminates the kinematic segrega- 
tion of galaxies with different ages often reported in litera- 
ture. These two classes of galaxies are clearly separated in 
the {g — i)-i color-magnitude diagram (Fig. 6), and most D99 
k-type galaxies belong to a well-defined cluster red sequence. 

3.1. Velocity Distribution 

For each cluster [On] emitter we computed its cluster- 
centric line-of-sight velocity (or simply radial velocity) in the 
rest-frame of cluster: 

. _ Z-Zc\ 

Avios = c- , 

1+Zcl 

where z is the observed redshift of the galaxy and za = 0.4069 
for Abell 851. Fig. 7 shows the velocity distribution of 
cluster [O ll] emitters and D99 cluster member galaxies. It 
is quite obvious that the velocity distribution of our clus- 
ter [On] emitters is not well described by a Gaussian. A 
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Fig. 6. — The (g-i)-i color-magnitude diagram of all MLFOO objects de- 
tected in their broadband images (gray dots), D99 k-type galaxies (open di- 
amonds), and cluster [Oil] emitters (symbols as in Fig. 4). The g — i colors 
and magnitudes of D99 objects are determined by those of spatially cross- 
correlated objects in the MLFOO catalog, which includes the g - i colors of all 
the objects detected in their broadband images. In the color version, the D99 
spectral types for the [O II] emitters (§4) from the best estimates of the [O II] 
and H(5 equivalent widths are indicated by blue [e(b)], green [e(c)], red [e(a)], 
yellow (k), and magenta (k-l-a/a-Hk). [See the electronic edition of the paper 
for a color version of this figure. \ 
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Fig. 7. — (Top) Cluster rest-frame line-of-sight velocity distribution of clus- 
ter [O II] emitters. The solid histogram corresponds to the spectroscopically- 
confirmed cluster members with [O II] emission, and the short-dashed his- 
togram is that of cluster [O II] emitters within the common survey area with 
D99. Top axis corresponds to the central wavelength of redshifted [Oil] 
emission line. The MLFOO on-band filter transmission curve is shown in a 
gray dash-dotted line with the transmission efficiency indicated on the right 
axis. (Bottom) The velocity distribution of the confirmed Abell 851 mem- 
bers with a k- (solid histogram) and k-Ha/a-l-k-type spectrum (long-dashed 
histogram) from D99. 



virialized population is expected to have a Gaussian veloc- 
ity distribution around the cluster systemic velocity. A K-S 
test on the distribution of [O ll] emitters rejects a Gaussian 
at a very high confidence level (> 99%). A peculiar fea- 
ture of the velocity distribution of cluster [Oil] emitters is 
the substantially small fraction in low-velocity [Oil] emit- 
ters (e.g., -500 km s"' < Avi^, < km s"'). Whether D99 
k-type galaxies show a similar deficit is not clear; statisti- 



cal tests involving the D99 sample are not very illuminating 
as its small sample size generally gives inconclusive results. 
For example, a two-sample K-S test of the velocity distribu- 
tions of all the [O II] emitters and D99 k-type galaxies with 
-2500 km s"' < Avi^, < -1-2500 km s"' gives - 92% con- 
fidence that the two distributions are drawn from the same 
parent population. Yet, if we use the subsamples restricted 
spatially within the common survey area of -3 .'74 < ARA < 
-1-3 f97 and -2f42 < ADEC < -l-l.'91,^ the same statistical test 
rejects at the 97% confidence level the hypothesis that they are 
drawn from the same parent distribution. At least D99 k-type 
galaxies do not show as pronounced a depletion near the clus- 
ter systemic velocity and might represent a virialized popula- 
tion; this is not a baseless assumption, since the D99 k-type 
galaxies appear to form a well-defined cluster red-sequence 
(Fig. 6), and such galaxies are likely to trace gravitationally- 
bound substructures (although it is almost certainly the case 
that Abell 85 1 has not been virialized over the entire cluster 
scale). 

The bi weight estimator and the jackknife method are used 
to estimate the velocity dispersions cr,, and their uncertainties 
(Beers et al. 1990): 1 160 ± 140 km s"' and 1010 ± 50 km s"' 
for D99 k-type galaxies and cluster [O ll] emitters, respec- 
tively. We must note, however, that Fig. 7 shows that a 
significant number of cluster [Oil] emitters with Avios < 
-1000 km s"' escaped detection by the MLFOO narrowband 
search due to the sharply dropping transmission efficiency. 
This would lead to a severe underestimate of the true veloc- 
ity dispersion for the cluster [O ll] emitting population; their 
velocity dispersion should therefore be taken as a lower limit. 

A simple kinematical treatment of infalling and virialized 
cluster galaxies in a cluster-scale potential leads to \T /V\ « 1 
and \T/V\ « 1/2, respectively, where T and V are the ki- 
netic and potential energies of galaxies. In this picture, the 
velocity dispersions of infalling and virialized galaxies are 
expected to follow the relation (T,,,infaii ~ 2'''^(T,,,vir. Having 
a higher velocity dispersion compared to that of virialized 
cluster members therefore is often invoked as a basic attribute 
of an infalling galaxy population. Previous studies generally 
found such kinematic segregations of cluster galaxies, with 
star-forming galaxy populations having higher velocity dis- 
persions (e.g., Colless & Dunn 1996; Dressier et al. 1999; 
Diaferio et al. 2001), although Rines et al. (2005) found Ut- 
ile evidence for the generality of this trend using their sample 
of emission-dominated and absorption-dominated galaxies in 
virialized regions of clusters. Such disagreements only illumi- 
nate the intrinsic difficulty of minimizing the effects of biases 
in cluster-scale surveys of galaxies, which require a substan- 
tial amount of spectroscopy. 

In short, the cluster [O ll] emitters appear not to be virial- 
ized, since their velocity distribution does not follow a Gaus- 
sian around the cluster systemic velocity. Although we can- 
not rigorously conclude due to the small sample size of D99 
galaxies, their velocity distributions restricted to the common 
survey area show some evidence of kinematic segregation be- 
tween the cluster [Oil] emitters and the D99 k-type galax- 
ies. The non-virialized state of [Oil] emitters naturally arises 
when they are at the infalling phase and have not spent a long 
time within the cluster gravitational potential. 

3.2. Spatial Distribution 

^ The coordinates are given with respect to 9''42'"56?2, 46°59'12"0 
(J2000). 
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Fig. 8. — Spatial distributions of cluster [O II] emitters (bottom row) and D99 galaxies {top row) with k- (open triangles) and k+a/a+k-type (filled triangles) 
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ically with respect to the cluster dynamical center Some of 
the spatially non-clustered [O ii] emitters may be field inter- 
lopers. 



Fig. 8, on the other hand, shows clear evidence for a kine- 
matic segregation in terms of spatial distributions of the clus- 
ter [O II] emitters and D99 k-type galaxies. In contrast to D99 
k-type galaxies, [O ll] emitters are only sparsely distributed in 
the dense, cluster core region; there are very few [O II] emit- 
ters inside or near the contour of the highest local galaxy den- 
sity. This trend holds true in aU the velocity bins in the figure. 
In Fig. 1 we have shown that the [O ii] emitter fraction is a 
strong function of local galaxy density, and Fig. 8 only visu- 
alizes the well-established empirical relation. In wide-field 
imaging of Abell 851, the red-sequence galaxies are shown 
to trace similar large-scale structures and substructures as ob- 
served via the photometric redshift technique (Kodama et al. 
2005), so it is of no surprise that the D99 k-type galaxies are 
found abundantly in dense regions. 

The very skewed distribution of blueshifted (Avios < 
km s"') [On] emitters is also evident in the southern 
(ARA ~ -K)' and ADEC ~ -6') and north-western regions 
(ARA ~ -A' and ADEC ~ +4'). Systematically redshifted 
[O II] emitters are present just in the north of the cluster cen- 
tral region as well. Abell 85 1 is known for the supposed fil- 
amentary structures extending out from its central region in 
these directions (KOI). The skewed distribution of Avios in- 
dicate the [On] emitters are part of infalling substructures, 
especially when spatially clustered as seen clearly in the fig- 
ure. 

In short, we find that the cluster [O ii] emitters tend to avoid 
the densest envirormients, but their skewed velocity and spa- 
tial distributions strongly suggest they are moving systemat- 



3.3. Infalling or Backsplash Population? 

True positions of the galaxies relative to the cluster cannot 
be determined from their redshift alone, however, and in ad- 
dition to the possibility of seeing field interlopers, even clus- 
ter members beyond the virial radius (/?vir) could also have 
gone through the cluster central region once or more in the 
past, constituting a so-called backsplash population (Balogh 
et al. 2000; Mamon et al. 2004). According to the numerical 
simulations, about half of the galaxies within 1 < R/Ryn < 2 
can belong to a backsplash population by the present epoch, 
z = 0. The backsplash scenario has been invoked to explain the 
presence of H i deficient galaxies beyond virial radii of nearby 
galaxy clusters, where the galaxies appear to be ram-stripped 
of their H I supplies during the passage through dense intra- 
cluster gas of cluster cores (Solanes et al. 2002; Sanchis et al. 
2002). Since the inferred physical mechanisms that transform 
cluster galaxies vary for infalling and backsplash populations, 
finding their dynamical state helps constraining SFHs of clus- 
ter galaxies. 

The simulations carried out by Gill et al. (2005) have shown 
that backsplash and infalhng populations exhibit kinemati- 
cally distinct velocity distributions beyond virial radii, where 
a backsplash population has about a factor of two smaller 
cluster-centric velocities than infalling satellites. They have 
also shown that the distribution of cluster-centric line-of-sight 
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Fig. 9. — (Top) Distribution of |Avios| for cluster [Oil] emitters (solid 
histogram) and D99 k-type galaxies (dashed histogram). The top axis shows 
the velocity normaUzed to an estimate of cluster velocity dispersion cTy = 
1260 km s~^ The fractional abundance is normalized to the total nvmiber of 
galaxies within each sample to faciUtate comparison. Lack of low-| Aviosl 
galaxies is a signature of infalling population (cf , Fig. 8 of Gill et al. 2005). 
We note that the samples are not restricted to 1 < i?/i?vir < 2. (Bottom) The 
velocity distributions of the cluster [O II] emitters inside (short-dashed) and 
outside (dash-dotted) the fiducial virial radius of 2 Mpc. The center of the 
cluster was taken to be at 9''43"03?3, 46°59'26;'5 (J2000), the brightest 
point in X-ray surface brightness from De FiUppis et al. (2003). 

velocities like Fig. 7 could imply the presence or absence of 
a backsplash population. Along this line, we show Fig. 9 to 
facilitate a comparison to Fig. 8 of Gill et al. (2005). The 
depletion of low- 1 Avios | [O ii] emitters shows quite a remark- 
able similarity to the distribution expected of an infalling pop- 
ulation in their simulation. It should also be noted that un- 
like Fig. 8 of Gill et al. (2005) the sample used in the top 
figure is not restricted to those within the projected radii of 
1 < R/Ryii < 2, where the difference would be more pro- 
nounced between backsplash and infalling populations. In the 
same figure, the D99 k-type galaxies also appear to show a 
peak at a higher velocity, the signature distribution of infalling 
population, just as the cluster [O ll] emitters do. However, we 
believe there are important differences, since the D99 k-type 
galaxies are drawn from a much smaller field of view (Fig. 8), 
are significantly brighter (Fig. 6), and are not necessarily rep- 
resentative in a statistical sense without applying proper cor- 
rections (Poggianti et al. 1999). We also lack D99 galaxies 
at large projected cluster radii. Furthermore, the peak around 
I AviosI — 1000 km s~' do not appear to indicate a real over- 
abundance of general galaxy population around that veloc- 
ity, since they are not spatially coincident; see Fig. 7 for the 
cluster central region sample. For these reasons, we believe 
the depletion of low-|Avios| galaxies is robust for the cluster 
[O II] emitters, yet the small sample size limits our ability to 
interpret the distribution of D99 k-type galaxies. 

In the bottom of Fig. 9, we show the velocity distributions 
of I O II] emitters inside and outside the fiducial cluster virial 
radius of 2 Mpc, the latter enabling a more direct comparison 
to the Gill et al. simulations. We caution that the distribu- 
tions shown here are rather sensitive to the choice of the fidu- 
cial virial radius, since the projected distribution of the cluster 
[O II] emitters are highly inhomogeneous as discussed earlier. 



Recently both Rines et al. (2005) and Pimbblet et al. (2006) 
carried out similar comparisons using composite galaxy sam- 
ples from nearby clusters, finding galaxies in infall regions are 
a mixture of infalling and backsplash populations. Since the 
Gill et al. results are for massive, moderately relaxed clusters 
at z ~ 0, these works produced more direct comparisons. In 
contrast to these authors' findings, the "peaks" of the veloc- 
ity distributions of our cluster [O ll] emitters extend toward 
higher velocities up to | Avios | /cr„ ~ 1 , whereas the galaxies in 
both Rines et al. and Pimbblet et al. show a steep decline in 
numbers at lAviosl/ci, > 0.5. For the R> 2 Mpc sample, we 
also observe a more pronounced depletion of low-vios ]Oll] 
emitters. These trends imply that the total [O ll]-emitting pop- 
ulation is probably a mixture of infalling and backsplash pop- 
ulations, yet the velocity distributions are generally broadened 
toward higher velocities than observed in local clusters, favor- 
ing a higher abundance of infalling objects. 

In short, although the survey area of our spectroscopic sam- 
ple is not as large as some of the most recent wide-field stud- 
ies of galaxy clusters, e.g., Treu et al. (2003) and Moran et al. 
(2005) go out to radii of ~ 5 Mpc in CL0024-I-1654, the anal- 
ysis of their kinematics suggests that a significant fraction of 
the cluster [O ii] emitters are likely to be still in the infalling 
phase. 

4. CLASSIFICATION BY SPECTRAL TYPE 

To probe transformation of galaxy properties, we need ob- 
servable indicators that are sensitive to the age of underlying 
stellar population. Classifying galaxies in terms of [O ii]-H5 
equivalent width plane is one way to put galaxies in the con- 
text of evolutionary sequence and has been used in literature 
(e.g.. Dressier et al. 1999; Balogh et al. 1999). The equivalent 
widths of [O II] emission and H6 absorption, respectively, are 
sensitive to the on-going star-forming activity and the pres- 
ence of intermediate-age stars, i.e., A- and early F-type star 
populations which live up to ~ 1 .5 Gyr with characteristically 
strong Balmer absorption lines. In the conventional interpre- 
tation of post-starburst galaxy spectra, the presence of strong 
Balmer absorption lines signals the recent cessation within 
the past ~ 1.5 Gyr of star-forming activity after which the 
intermediate-age stellar population makes a dominant contri- 
bution to the integrated optical/UV galaxy spectrum. There- 
fore, the [O II]-H5 plane and its variants have been useful for 
classifying a galaxy spectrum in relation to some mechanisms 
which halt or significantly reduce the amount of star-forming 
activity within a galaxy in its recent past. 

4.1. D99 Spectral Class Definitions 

We adopt the D99 spectral classification scheme to study 
the SFHs of cluster [Oil] emitters (Table 1).** We sfi-ess that 
by convention the measured H6 equivalent width in the [O ll]- 
Hd plane is the sum of emission and absorption equivalent 
widths. This means that Wo(H(5) < A does not mean the 
line is totally seen in emission; it merely means the strength 
of emission equivalent width is greater than that of absorp- 
tion at 116. Although we could separate the emission and ab- 
sorption components at US in many cases, we followed the 
conventional index definition to faciUtate comparisons to lit- 
erature. We used the equivalent widths measured using the 
flux-summing method for spectral indices (Paper I). We now 
briefly review the conventional interpretation of each spectral 

^ For a few objects the [O ll]A3727 or Hi5 line moved out of observed 
spectral range, in which case their spectra were classified into "N/A"-type. 
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TABLE 1 
D99 Spectral Classification 





Definition 


Cluster [O II] 


Emitters 




Wo([Oll]) 




N{S/N > 2f 


A'(all)'' 


Spectral Class 


(A) 


(A) 






e(b) 


<^0 


<+4 


3 


23 


e(c) 


^0 — 5 


<+4 


14 


29 


e(a) 


<-5 


>+4 


40 


53 


k+a/a+k 


>-5 


>+3 


2 


4 


k 


>-5 


<+3 


1 


2 



Note. — The spectral indices Wo([0 ll]) and Wo(H5) are mea- 
sured with the flux-summing method (Paper 1). 

" The number of cluster [O II] emitters for which the classification 
is done with the equivalent width measurements at S/N > 2. 

^ The number of cluster [Oil] emitters with each classification 
from the fuU sample. 

class. For more detail, readers are referred to the original pa- 
pers by D99 and Poggianti et al. (1999). 

e(b)-type. — Galaxies in this class have a strong [O ll] emis- 
sion, indicative of vigorous on-going star-forming activity. 
Low-lying recombination lines of the Balmer series may also 
be seen in emission. The duration of a starburst cannot be 
longer than about a few 100 Myr, however, since after that 
timescale the buildup of continuum flux near the line de- 
creases the [O ll] emission equivalent width to a degree that 
the spectrum would be classified into e(c)-type. 

e(c)-type. — Typical e(c) galaxies are a star-forming galaxy 
with only a moderate star-forming activity at a relatively con- 
stant rate. Local late-type spirals typically fit in this category. 

e(a)-type. — Typical galaxies in this class have strong 
Balmer absorption lines as well as moderate to substantial 
[Oll]A3727 emission. Based on the comparison with a lo- 
cal sample and high degree of dust extinction implied by low 
[O II] /Ha line ratios, Poggianti & Wu (2000) have associated 
typical e(a)-type galaxies to dusty starburst galaxies in which 
the light from the youngest starbursting population is much 
more extincted than that from the intermediate to old age stel- 
lar populations (see also Poggianti et al. 2001). 

k+a/a+k-type. — So-called post-starburst galaxies (also 
called E+A or H(5-strong galaxies) with little or no [O ll] emis- 
sion and strong Balmer absorption lines belong to this class. 
This class of galaxies is often studied as a missing hnk be- 
tween passively-evolving galaxies in clusters and actively- 
star-forming galaxies abundantly observed outside clusters. 

k-type. — Passively evolving, non-star-forming galaxies 
like local ellipticals belong to this class. A typical k-type 
spectrum is dominated by the Ught from K-type star popu- 
lation, hence the designation. 

A few representative spectra, actually observed in our spec- 
troscopic survey, are shown in Fig. 10 with their correspond- 
ing D99 spectral class. Table 1 also lists the numbers of clus- 
ter [0 11] emitters in each spectral class. As expected, most 
[O II] emitters belong to one of the emission-line classes [i.e., 
e(b), e(a), and e(c)], but six objects, all of which do show 
weak [On] emission, are classified into non-emission-line 
classes. By making a cut at Wo([0 11]) > -5 A, the D99 classi- 
fication scheme does not strictly remove the contamination of 
emission-hne objects into their "passive" classes. Among the 
four k-na/a-nk type objects, three are strong Balmer-absorption 
systems with residual [O ll] emission. One of the two k-type 
objects appears to have an early-type galaxy spectrum but also 
has [O II] and Balmer emission fines. The remaining two ob- 
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Fig. 10. — A set of observed high-S/N spectra classified into each D99 
spectral class. The regions affected by very poor sky subtractions are masked 
from each spectrum in this figure. Note that a bright skyline and telluric 
absorption lines generally contaminated redshifted H/3 absorption continua, 
substantially increasing the uncertainties in the H/3 line flux measurements. 

jects appear genuine ambiguous classifications due to poor- 
S/N spectra. 

4.2. Distribution of [On] Emitters by Spectral Class 
For our sample of cluster [O ll] emitters, the limiting rest- 
frame [Oil] equivalent width is ~ -8 A at i ~ 23, so there 
are few galaxies of k- or kn-a/aH-k-type, which must have 
Wo([Oll]) > -5 A by definition; most [O ll] -selected galax- 
ies belong to one of the emission-line galaxy types: e(b)- 
, e(c)-, and e(a)-types. In Fig. 1 1 we show the distribution 
in the [Oii]-H(5 plane for Abell 851 cluster [On] emitters 
with Wb([Oii]) < -8 A and i<23. To see how the disfii- 
bution compares to "field" [O ll] emitters, we measured the 
same spectral indices using the fourth data release of the 
Sloan Digital Sky Survey (SDSS; Adelman-McCarthy et al. 
2006), the Canada-France Redshift Survey (CFRS; Lilly et 
al. 1995; Hammer et al. 1997), and data release one of the 
DEEP2 Galaxy Redshift Survey.^ For each survey, we only 
included galaxies with Wb([0 ii]) < -8 A to compute the rel- 
ative abundances of spectral types. Table 2 summarizes the 
basic characteristics of the samples. The typical size of the 
cluster [O ll] emitters is 2"-3" (~ 13 kpc at z = 0.4), observed 
with 1 "-width slit in our spectroscopic program. The spectro- 
scopic apertures in SDSS (3" diameter fiber), CFRS (l."75- 
width sUt), and DEEP2 (1 "-width slit) sample the regions of 
galaxies roughly within 5-6 kpc, > 1 1 kpc, and > 8 kpc, re- 
spectively, at the median redshift of each survey. Although 
we do not make an attempt to correct for the differences, the 
regions sampled by the surveys are comparable or at least 
greater than our spectroscopic aperture. All the surveys are 
magnitude-limited to different depths, so field samples can 
only be constructed to a depth brighter than the completeness 
limit of the cluster [O ll]-emitting sample. 

We constructed rough volume-limited samples by scaling 
the number of galaxies within similar rest-frame B-band ab- 
solute magnitudes by suitable luminosity functions from ex- 
isting Uterature. We used software package kcorrect (v4_l_4; 
Blanton et al. 2003) to obtain absolute magnitudes from a 

' http://deep.berkeley.edu/ 
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Fig. 11. — The distributions of [O II] and Hd equivalent width for Wb([0 II]) < -8 A objects for magnitude-limited samples. See § 4.2 for the selection methods 
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and their size is proportional to E(B-V). The \a error bars indicate the median value for a measurement within each sample. 



TABLE 2 

Comparison Samples of Field [On] Emitters 



Survey z Range Magnitude Cut Luminosity Function Number of Spectra 

SDSS(DR4)'' 0.09-0.11 i < 23'' Ball et al. (2005)'^ 41613 

CFRS'' 0.2-0.7 /(ABX22.5 N/A'= 51 

DEEP2(DRl/ 0.8-1.0 R(AB)<24.1 Willmer et al. (2005)8 1492 

Note. — The spectral indices Wo([0 ll]) and Wo(H5) are measured with the flux-summing 
method (Paper I). Each field sample consists of galaxies with Wo([0 II]) < —8 A. 

" Only extended sources which are spectroscopically classified as galaxy with zconf > 0.85 
are included. 

The magnitude limit of SDSS survey is / Ri 21.3; we include fainter galaxies to increase 
the number of galaxies near Mb » -19. All photometry in our SDSS samples are Petrosian 

magnitudes. 

eClass PCA bivariate luminosity function (their Fig. 5). Within each magnitude bin and 
eClass range, an additional correction for the fraction of galaxies with Wo([0 II]) < -8 A is 
applied to remove the contribution to the luminosity function from non-emission-line galaxies. 

Only objects with their confidence class 0^ (Le Fevre et al. 1995) are included. As the 
CFRS spectra were obtained at considerably different spectral resolutions (40 A) than other sur- 
veys (~ 3 A for SDSS and ~ 1 A for DEEP2), we used slightly different windows of 3653 A- 
3706 A, 3706 A-3748 A, and 3748 A-3801 A ([Oil]) and 4030 A-^078 A, 4078 A^126 A, 
and 4126 A— 4170 A (H5) for the blue continuum, fine, and red continuum, respectively. 

° No attempts have been made to correct the sample using a luminosity function. 

^ Only objects with zquality > 3 are included. 

8 Schechter function for blue galaxies (their Table 4). 



set of available broadband photometry (corrected only for 
Galactic extinction) for each galaxy survey. The cluster 
[O ll] emitters near the MLFOO narrowband survey limit have 
MB(Vega) ~ -17.5 (before correcting for internal extinction). 
Since the completeness severely decreases toward that magni- 
tude, we limited our comparisons roughly to the survey limit 
of our DEEP2 sample, Mb < -19; for both the SDSS and 
DEEP2 samples, the rough completeness limit in terms of ab- 
solute magnitude is Mb ~ -21, so the uncertainties in galaxy 



fractions increase significantly toward the faint end of lumi- 
nosity function for that reason. We caution that the above 
procedure was not followed for the CFRS sample because 
of limited availability of photometric data and significantly 
lower-resolution spectra; thus the CFRS data points should 
only be used for a rough comparison. We further note that 
using purely magnitude-limited samples did not significantly 
change the overall trends to be discussed now. 
Looking at the SDSS field sample in the [O ii]-H(5 plane 
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Fig. 12. — Fractional abundance of each D99 spectral class relative to all 
Wo([0 II]) < -8 A galaxies in each survey. FUled symbols denote the Abell 
851 [O II] emitters. Each sample is shown roughly at its median redshift with 
a slight offset for clarity. 

is illuminating because of its large sample size and relatively 
small measurement uncertainties. We notice the distribution 
is dominated by typical star-forming galaxies [e(c)] , with rela- 
tively small pure starbursting population [e(b)]. The galaxies 
with e(a)-type spectrum appear to form some extension, to- 
ward stronger HS absorptions, of general star-forming galaxy 
population with e(c)-type spectrum. The DEEP2 field sam- 
ple generally suffers from large uncertainties; the distribution 
of galaxies with low [O ii] equivalent width does not appear 
to extend well down below ^ -10 A in contrast to the SDSS 
sample, which indicates more vigorous star-forming activity 
expected of galaxies at high redshift. 

In Fig. 12 we show the relative abundance of each D99 
class within each sample of Wo([Oll]) < -8 A galaxies at 
different redshift. As expected, the e(b)-type galaxy frac- 
tion was higher at high redshift, reflecting the fact that trig- 
gering of young starbursts was more frequent as the universe 
had a greater amount of star-forming activity out to z 2. 
A similar redshift trend is seen in the e(a)-type fractions in 
the field galaxy samples, but in Abell 851 the e(a) fraction is 
higher than or at least comparable to that in the field within 
the probed redshift range. It is important to note that the plot 
shows the fractional abundance of different spectral classes 
within each W()(LOll]) < -8 A galaxy samples. Since the 
[O II] emitters constitute a progressively smaller fraction of 
galaxies toward lower redshift (Nakata et al. 2005), the ob- 
served trend of e(b) and e(a) galaxies in the field is stronger, if 
the fractions are computed within the full galaxy samples. For 
the same reason, the e(c) fractions within the full samples do 
not evolve very much in the field, if we adopt the [O ll] galaxy 
fraction in Nakata et al. As noted earher, e(a)-type spectrum 
has been attributed to dusty starburst galaxies. If this identifi- 
cation is true, even the e(a) galaxies constitute a starbursting 
population, implying that starburst phenomena in Abell 851 
are at least as frequent as z ^ 1 field [O ll] emitters. 

In computing fractional abundance of each emission-hne 
galaxy class, we opted to take each data point in the [O ii]-H^ 
plane as a two-dimensional Gaussian probability distribution 
with its best estimates and standard deviations given by the 
measurements of [O ii] and H6 equivalent widths and their 
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Fig. 13. — The [On] and H<5 equivalent widths as a function of age for 
constant burst (bottom) and exponentially-declining SFH (top) models with 
the fiducial parameters E(B-V) = 0.5 and £*/£g = 0.44 (Appendix A). The 
models were run for star-forming timescales of r = 10 Myr (solid line with 
square symbols), 100 Myr (dashed with circles), 1 Gyr (short-dashed with 
triangles), and 2 Gyr (dotted with diamonds) models. The measurements 
are made at the ages (1, 2, 5) X (10^, 10*, lO' yr) and are indicated by 
symbols. The dashed gray line indicates the boundary of burst and non-burst 
populations, i.e., e(b) and e(c), as defined in terms of [O II] equivalent width. 
The dash-dotted grey line indicates the threshold [Wo(H<5) > 4 A] for e(a) 
population. 

uncertainties. Each data point integrates to unity over the en- 
tire [Oll]-H(5 plane, but low-quality measurements give rise 
to some probability for each point to leak into other spectral 
classes than that designated solely by their best estimates. A 
significant number of both AbeU 85 1 and field [O ii] emit- 
ters have spectra with low S/N — about a few per pixel for 
Abell 851 (Paper 1) and even smaller for the DEEP2 sam- 
ple; by following the above procedure, we avoid giving too 
much confidence to the classification using the best estimates 
of such low-S/N measurements. Dividing the sum of inte- 
grated probabilities in each spectral class by the total number 
of galaxies in the [Oll]-H5 plane yields the fractional abun- 
dance. For Abell 85 1 sample, we have taken into account the 
spectroscopic target selection bias by multiplying each g - i 
color subsample by an appropriate correction factor (Paper I). 

5. STAR FORMATION HISTORY OF CLUSTER [O II] EMITTERS 

Analyzing the cluster [O ll] emitters in terms of possible 
evolutionary scenarios requires modeling their spectral evo- 
lution using stellar population synthesis models. The D99 
spectral classification scheme employs a combination of stel- 
lar and nebular lines that are sensitive to the time evolution 
of galaxy spectra over a wide range in stellar age. Hence 
our spectral modeling considers the time evolution of syn- 
thetic spectra of underlying stellar populations as well as cor- 
responding nebular emission. Furthermore, combining these 
two components enabled us to incorporate the effect of dif- 
ferential dust extinction on stellar and nebular components, 
which appears essential for reproducing e(a) spectrum (Pog- 
gianti et al. 2001). The detail of our modeUng scheme is out- 
lined in Appendix A. 

5.1. Time Evolution in the [Oll]-H5 Plane 

The [O II] emitters with e(a)-type spectrum make up a sig- 
nificant fraction of the cluster [O ll] emitters (Fig. 1 1), so we 
place an emphasis on reproducing the e(a)-type spectrum for 
its dominance in our sample. We immediately excluded in- 
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stantaneous burst models, which were also modeled, as they 
generally fail to reproduce a typical e(a) spectrum with strong 
H(5 absorption and [O ii] emission at the same time. The val- 
ues E{B-V) = 0.5 and Et:/Eg = 0.44 are chosen as fiducial 
numbers representing our [O ll] emitting sample. 

In Fig. 13 we show the evolution of [On] and H(5 equiva- 
lent widths of constant burst and exponentially-declining SFH 
models. The 116 equivalent width cut for e(a)-type galaxies is 
indicated by the horizontal dash-dotted lines. Strong H5 ab- 
sorption is generally achieved by the models with a shorter 
star-forming timescale, t < 100 Myr, and the phase can con- 
tinue up to a few Gyr, roughly comparable to the lifetime of 
intermediate age stellar population. These short-r galaxies 
are generally classified into e(b)-type for up to ~ 100 Myr 
during which they show strong [Oil] emission (< -40 A), 
while their HS absorption has not had a time to grow as much. 
After the vigorous star formation, the short-r galaxies go into 
a post-starburst (k-na/a-nk) phase rather quickly due to their 
fast fading [O ii] emission. Therefore short-r models with 
the fiducial parameters of E(B-V) and E^/Eg are unlikely 
to explain the abundance of e(a)-type galaxies in our sample 
(Fig. 11), since the short-r models stay in the e(a) phase for 
no more than ~ 100 Myr. 

Strong [O II] emission and US absorption equivalent widths 
are achieved for a more extended duration by the models with 
star-forming timescales r ~ 1 Gyr. A major difference be- 
tween the constant burst and exponentially-declining SFHs 
with similar r is how star formation continues after f ~ r 
(which characterizes roughly the "burst" duration). Since 
O- and B-type stars can totally dominate the light from a 
galaxy when present, the H6 absorption line grows substan- 
tially only after a significant amount of O- and B-type stars 
are gone. Since star formation totally ceases after a burst, the 
[O II] emission equivalent width from constant burst models 
quickly goes below the boundary for emission-line galaxies 
[Wo(|0 III) > -5 A] after a few Myr. On the other hand, the 
r ^ 1 Gyr models with an exponentially-declining SFH can 
spend up to several Gyr within the emission-line categories, 
which is desirable for reproducing a high abundance of e(a) 
galaxies. 

Now, Fig. 14 illuminates the effect of dust extinction in 
the [Oll]-H(5 plane with a set of r 1 Gyr exponentially- 
decUning SFH models. A range of reddening E(B-V) = 0.0- 
1.5 roughly spans the spread of [O ii] emitters. We point out 
that the equivalent width is a measurement insensitive to the 
presence of dust reddening when the continuum and line flux 
are both attenuated by the same amount, in which case the 
model tracks do not deviate from the dust-free, E(B-V) = 
model. Therefore we introduced the effect of differential ex- 
tinction by varying E^/Eg, the ratio of E{B-V) in steUar to 
nebular component. It is clear that differential extinction is 
essential for pushing a model track into the e(a) regime, as 
the H(5 equivalent width for the dust-free track turns over at 
~ 4 A, barely crossing the boundary for the e(a) phase; tak- 
ing E^/Eg < 1 extincts the nebular HS emission more than the 
stellar continuum around it, enabling the equivalent width to 
grow slightly larger in absorption (i.e., less emission filling). 

However, the [On] emission, being a bluer line, suffers 
even more significantly from dust extinction as well. This 
in turn leads to a degeneracy between the age and the amount 
of reddening in the [Olll-H(5 plane. For example, a highly 
reddened E{B-V) =1.5 starbursting galaxy at the age of 
100 Myr is hard to distinguish from a E(B-V) = 0.5 galaxy 



TABLE 3 

i'-logS SUBSAMPLES 



Category logS N 



HB 


> 1.6 


<21 


10 


TB 


0.0-1.6 


<21 


3 


LB 


<0.0 


<21 


6 


HVI 


> 1.6 


21-23 


14 


TM 


0.0-1.6 


21-23 


29 


LM 


< 0.0 


21-23 


6 


HF 


> 1.6 


> 23 


6 


TF 


0.0-1.6 


>23 


19 


LP 


<0.0 


>23 
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at 3 Gyr, since they both would have Wo([Oll]) ~ -11 A 
and W()(HJ) ~ 5 A. The inferred age of star-forming activity 
can be quite different without the knowledge of reddening for 
breaking such a degeneracy. 

5.2. Star Formation History by Composite Spectra 

We recall that it was largely the interpretation of low-S/N 
HS equivalent width measurements which led to the contro- 
versy regarding the frequency of starbursts in cluster galaxies 
(Dressier et al. 1999; Balogh et al. 1999). Although our anal- 
ysis of e(a)-type galaxies in the [Oii]-H(5 plane in the pre- 
vious section independently confirms the high abundance of 
cluster galaxies with strong Balmer absorption fines, also re- 
ported recently by Dressier et al. (2004), the relatively low 
spectral S/N of faint [0 11] emitters limits our ability to con- 
striiin the individual SFHs with much confidence. For our typ- 
ical spectra, a 3a detection of HS absorption equivalent width 
at Wb(Hi5) ~ 5 A, for example, needs per pixel S/N ~ 10 in 
the continuum, a requirement not quite met by most of our 
spectra (Sato & Martin 2006). We therefore supplement our 
analysis by composite spectra, which smear out the SFHs of 
individual galaxies but would help us find the characteristic 
SFHs of subsamples. 

Our scheme for subsampling cluster [O ll] emitters is moti- 
vated by the observational evidence that the transformation of 
galaxies depends both on their intrinsic property and on their 
local envirormient. We chose i magnitudes and local galaxy 
densities E to characterize these attributes; Fig. 5 shows the 
classification scheme, and Table 3 shows the number of ob- 
jects within each subsample. We effectively make subsam- 
ples residing in the local environments characteristic of low- 
density, field-like regions (L; logS < 0.0), the transitional 
zone (T; 0.0 < logS < 1.6), and the high-density regions (H; 
logS > 1.6). The highest density cut was made close to the 
KOI threshold density. Each density bin is further divided 
into three / magnitude bins: bright (B; / < 21), medium (M; 
21 < J < 23), and faint (F; i > 23). The subsampling yields 
composite spectra of quafity ranging from per pixel contin- 
uum S/N ^ 3 for faint subsamples to S/N ^ 30 for bright sub- 
samples. For faint subsample this is still an improvement, as 
most ;■ > 23 [O ll] emitters have S/N < 2. 

The goal of this exercise is to study the SFHs inferred from 
representative composite spectra of various cluster [Oii]- 
emitting subsamples by placing them in the [Oll]-H(5 plane. 
For this reason the method for coadding spectra must not ar- 
tificially change the resulting |0 ll| and Hi5 equivalent widths 
of composite spectrum, i.e., continuum flux strength relative 
to the fine flux must be scaled properly at these lines. Averag- 
ing spectra normalized to a common mean flux density over a 
relatively featureless continuum region for the entire spectrum 
is not desirable, especiaUy for [O ii]A3727, since the emission 
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Fig. 14. — The effect of dust extinction on the observed spectral indices of a fiducial model track (exponentially-declining SFH with r = 1 Gyr, E(B—V) = 0.5, 
and Et,/Ej, = 0.44) indicated by circles connected by dashed line. (Left) The output spectra from the model were extincted by a different amount of E(B-V) at 
constant £» /Eg = 0.44. (Right) Stellar and nebular continua were differentially extincted at a ratio given by £« /Eg at constant E(B-V) = 0.5. Filled data points 
indicate models at 10 Myr, 100 Gyr, and 1 Gyr. Filled gray circles are observed cluster [O II] emitters. 

line appears within a sloping continuum region of the Balmer 
break. Instead we normaUzed spectra to the continuum flux 
density at line center of each line of interest, after the spec- 
tra are shifted to their rest frame wavelengths. All spectra are 
rebinned to a common dispersion at 1 A per pixel before av- 
eraging to produce the final composite spectrum. A potential 
bias may be introduced to a composite spectrum depending 
on how each spectrum is weighted or not weighted against 
others within the common subsample. We cannot weight by 
such quantities as inverse variance, since it is strongly depen- 
dent on source luminosity, an intrinsic property of subsam- 
ple population which we wish to control. In order to treat 
each [O ii] emitter within a subsample on equal footing, the 
composite spectrum was computed without weighting. The 
results are also checked against the composite spectra made 
from medianing, which is generally more robust against out- 
liers. We also excluded objects identified to be AGN-like from 
our composite samples. 

In Fig. 15 we see where the subsamples lie in the [Oll]- 
H(5 plane, overlayed on the evolutionary tracks of our fidu- 
cial models for E{B-V) = 0.0- 1.5. Within similar local en- 
vironments, the [O II] equivalent width of [O ii] emitters is 
strongly dependent on their i magnitude; luminous galaxies 
tend to have a weaker [O ll] emission, in line with what we 
found in Paper 1. Remarkably, their H6 equivalent widths are 



-10 



e(b) 


e(a) 

^ ' - 




O LF / 
LM / 

♦ LB /■-, 
O TF / * 
•a TM / 

• TB / 

- ^ HF efcl / / 

" SB , 1 ,/ , ,/ 


^/ > 



w^m (A) 



Fig. 15. — The [O II] and H5 equivalent widths from the composite spec- 
tra of !-S subsamples. The Icr uncertainties were computed formally from 
the sigma vectors of spectra to be coadded. The gray lines indicate our fidu- 
cial emission galaxy models at E(Ji-V) - 0.0 (solid), 0.5 (long-dashed), 1.0 
(short-dashed), and 1.5 (dotted); see Fig. 14. The uncertainties in Wb([Oll]) 
for some data points are smaller than their symbol size. 
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Fig. 16. — [O II] and HS equivalent widths of the composite spectra from 
the i-E(B-V) subsamples (§ 5.2). A cut was made at E(B-V) = 0.5 to divide 
[On] emitters into galaxies with low and high reddening, and each subsample 
was further divided into three luminosity bins with cuts at / = 21 and 23. The 
uncertainties in Wo([0 II]) for some data points are smaller than their symbol 
size. 

nearly constant at 4-5 A, and only the LF sample stands out 
with a small US equivalent width in absorption (~ A). The 
LF sample is classified into e(b)-type, yet its HS measure- 
ment has a large uncertainty; within 2<t, the LF data point is 
consistent with being e(a)-type. While this is true, it is gen- 
erally difficult for a spectrum to have strong |0 II] emission 
and absorption equivalent widths at the same time because 
of emission filling at HS. A very strong [O ii] emitter with 
Wb([0 ll]) < —100 A would always have HS equivalent width 
in emission, i.e., Wb(H(5) < A (Fig. 14). Therefore, with a 
reasonable confidence in the very strong [O ii] emission, it is 
hkely that the LF composite galaxy is a true e(b)-type galaxy. 

It appears our model tracks with E(B-V) = 0.0- 1 .5 almost 
entirely cover the range of measurements from the compos- 
ite spectra. While this may indeed imply the importance of 
dust extinction, we expect the i-T, sampling scheme to have 
degeneracy on the amount of dust reddening in the [O ll]-H5 
plane. Composite sampling makes more difficult to interpret 
the distribution seen in Fig. 15 in terms of the amount of red- 
dening, since we have seen there is a positive trend of E{B-V) 
on galaxy luminosity (Paper I). 

To better interpret the effect of reddening, we make an- 
other subsampling based on E(B-V) in Fig. 16. On a model 
evolutionary track, increasing E(B-V) has an effect of push- 
ing the entire track toward weaker [O ii]-emitting e(a)-type 
galaxy (Fig. 14). If the dust reddening is the primary param- 
eter for the spread of [O ll] emitters, resampling them based 
on E(B-V) might form a one -parameter family of data points 
along the evolutionary track on the [O Il]-H(5 plane, but we see 
this is not the case. A cut was made at E(B—V^ — 0.5 to divide 
[O II] emitters into galaxies with low and high reddening, and 
each subsample was further divided into three luminosity bins 
with cuts at ; = 21 and 23. For the subsamples with similar lu- 
minosities, the effect of reddening does not appear to be as 
strong as expected from the shift in the model tracks arising 
from varying the amount of reddening alone, except perhaps 
for the brightest subsamples. The general trend of decreas- 
ing strength of [O ll] equivalent width with luminosity is still 
observed. 
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E{B-V) 

Fig. 17. — Rest-frame [O II] equivalent widths of cluster [O II] emitters as a 
function of their reddening iiCB-V). Data points are divided into low- (i > 23; 
open triangles), intermediate- (23 > i > 21; half-filled circles), and high- 
luminosity (i < 21; fiUed squares) sources. Objects identified to be AGN-hke 
are excluded from the figure. 

The luminosity trend is particularly hard to interpret for 
moderate [O ii] emitters in the presence of dust, since evo- 
lutionary tracks give rise to degeneracy in burst age and the 
amount of reddening. For example, in Fig. 16 the luminous 
(B), E(B-V) > 0.5 subsample lies close to the point at which 
the model tracks with E(B-V) = 0.5 and 1.5 intersect. As 
far as the amount of reddening is concerned, both the evolu- 
tionary tracks are very roughly consistent with the measured 
E(B-V), yet the inferred burst age is very different depending 

on what the actual E(B-V) is 100 Myr for E(B-V) = 1 .5 

and ^ 3 Gyr for E(B-V) = 0.5. Rather than constructing com- 
posite samples, treating galaxies individually with distinct, 
more precise reddening measurements would better constrain 
their SFHs in this case, if the quality of spectra allows it. 

To better elucidate the coupled effect of luminosity and red- 
dening in the [Oll]-H(5 plane, in Fig. 17 we relate the [Oil] 
equivalent width, luminosity, and the amount of reddening. 
We clearly see that luminous (i < 21), weak-[Oii] emitters 
[Wb([0 II]) - -10 A] are highly reddened with E(B-V) - 1. 
The bright subsamples in Fig. 15, LB, TB, and HB, lie in 
the region roughly bounded both by the tracks at the age of 
~ 1-2 Gyr of low reddening E(B-V) ~ 0.0-0.5 model galaxy 
and the tracks at the age of ~ 100 Myr of high reddening 
E(B-V) ~ 1.0-1.5 model galaxy. The observed high redden- 
ing suggests the younger, highly reddened burst is a more nat- 
ural interpretation of luminous, weak-[Oll] emitters. In gen- 
eral, we see that large amount of reddening can make young 
starburst galaxies appear in the same loci as less reddened 
starburst at older burst age, and we tend to observe them in 
the e(a)-type spectrum. A high amount of reddening typically 
observed in the cluster |Oll| emitters \E{B-V) > 0.5] sug- 
gests that younger burst (< 1 Gyr) is more likely to explain 
the dusty [O ii]-emitting population. 

6. DISCUSSION 

Recent wide-field surveys of cluster galaxies have shown 
convincing evidence that the transformation of galaxy proper- 
ties starts well beyond the viriaUzed part of galaxy clusters, 
signaling the importance of local environments for prepro- 
cessing of cluster galaxies (e.g., Balogh et al. 1997; Treu et 
al. 2003; Rines et al. 2005; Moran et al. 2005; Pimbblet et al. 
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2006). The star-forming properties of Abell 851 [On] emit- 
ters have shown a strong dependence on local environment as 
often reported in Uterature; the [O ii] emitter fraction depends 
strongly on local envirormients (e.g., Balogh et al. 2004b; Ko- 
dama et al. 2004; Nakata et al. 2005), while the strength of 
[O II] equivalent width is relatively independent of environ- 
ment (Balogh et al. 2004a,b; Kodama et al. 2004; Rines et al. 
2005). These properties in [O ll] emission are consistent with 
the interpretation that the mechanisms that give rise to the 
suppression of star-forming activity either act on very short 
timescales or are sufficiently rare, i.e., the abundance of tran- 
sient objects is relatively small. Furthermore, there appears 
a luminosity segregation of the [O II] -emitting populations in 
dense environments — at the luminous end, there is an indica- 
tion that AGN activities could play a role in suppressing star- 
forming activity, whereas cluster [O ii] emitters are notably 
absent at the faint end. 

6. 1 . Timescales of Star-Forming Activity 

We put some constraints on the timescales of star-forming 
activity in the cluster [On] emitters using a combination 
of stellar population synthesis and nebular photoionization 
codes. Figures 14 & 15 show that the spread of cluster [O ll]- 
emitters in the [Oll]-H(5 plane may be covered well by 
exponentially-decUning SFH models of r ~ 1 Gyr, with vary- 
ing amount of reddening up to E{B—V) ~ 1.5. A wide range 
of reddening, however, gives rise to a degeneracy between the 
burst age and reddening — a young, very dusty starburst and a 
less dusty star-forming galaxy at an older stage of evolution 
can appear at a similar locus in the [O ii]-H(5 plane. 

The composite spectra of the subsamples in the j'-E plane 
indicate that the galaxy luminosity is a major factor contribut- 
ing to the spread of [Oil] equivalent widths (Fig. 15). This 
is expected, since there exists an empirical anti-correlation 
between the [On] equivalent width and galaxy luminosity 
(Paper I). The interpretation of equivalent width is somewhat 
complicated by the fact that it also is a measure of relative 
fraction in young and old stellar populations; a galaxy with a 
small vigorously star-forming regions spread within a much 
larger region of old stellar population would have a weaker 
[O II] emission equivalent width, for example. Our spectral 
modeling assumed simplistic SFHs for which 100% of galaxy 
mass is produced in a single episode of star formation. A more 
reasonable scenario involves patches of young star-forming 
regions embedded within older stellar populations. In such 
a case, both [On] emission and absorption equivalent 
widths from spatially integrated spectra would shift slightly 
toward the origin (W = A) in the [O ll]-IW plane, in which 
the amount of shift is sensitive to the relative abundances of 
young and old stellar populations. Albeit modeling extreme 
cases, it is still encouraging that a simple treatment of dif- 
ferential extinction, as suggested by Poggianti et al. (2001), 
have reproduced the spread of observed [O II] emitters in the 
[Oll]-H(5 plane. 

Fig. 1 1 shows among the most abundant [O ll] emitters in 
Abell 851 and in the field at high-z are those with the e(a)- 
type spectrum. We re-emphasize that properly taking account 
for the differential dust extinction in the stellar and nebular 
components of a galaxy spectrum is essential for interpret- 
ing the e(a)-type spectrum. The dust-free model track barely 
reaches the absorption equivalent width strong enough to 
be considered an e(a)-type in Fig. 14. The identity of H5- 
strong galaxies with moderate-to-strong [On] emission has 
been somewhat enigmatic, since emission-filling at ab- 



sorption line can reduce its apparent strength seen in absorp- 
tion. Although the nature of e(a) galaxies and its role in 
cluster galaxy evolution remain debatable, the original inter- 
pretation of e(a) galaxies as dusty starbursts by Poggianti et 
al. (1999) has recently been given considerable boosts. A 
study of SDSS galaxies at z ~ 0.1 by Balogh et al. (2005) 
has concluded that they are normal or dusty spiral galaxies 
with E(B-V) ^ 0.7, assuming the Milky Way extinction of 
Ry = 3.1; a recent study of mid-infrared sources also have 
supported the picture that e(a) galaxies are dusty starbursting 
systems (La Franca et al. 2004). 

Although not a direct proxy for the presence of A-type 
stellar population, a very large positive value of HS equiv- 
alent width must indicate a dominant A-type star contribu- 
tion over optical spectrum, with the US emission-filling po- 
tentially masking their presence by reducing the equivalent 
width seen in absorption. H(^-strong galaxies have been stud- 
ied extensively in the context of post-starburst/star-forming 
systems, since a substantial amount of young O- and B-type 
stars have to be gone before the intermediate-age stars can 
dominate the optical spectrum. The simplest SFH that pro- 
duces a post-starburst/star-forming spectrum is a truncation 
SFH after an episode of starburst or star formation. For an 
e(a)-type spectrum, however, some residual star-forming ac- 
tivity after a burst is necessary to reproduce [O ll] emission in 
dust-free models (Poggianti et al. 1999). Ourmodehngof e(a) 
spectra also does not favor simple truncation of star forming 
activity (§5.1). 

When differential dust reddening is accounted for, r ^ 
1 Gyr models reproduce well the spread of cluster [O ii] emit- 
ters in the [O ii]-H5 plane. That the timescale of r ~ 1 Gyr 
is of order of A-type star lifetime is not a coincidence; large 
116 absorption equivalent widths simply require A- or early- 
F-type star populations. This brief window of order 1 Gyr 
during which HS absorption is strong identifies either the re- 
cent cessation of starburst or on-going starburst. The former is 
the case for a post-starburst galaxy since the nebular emission 
fades quickly due to the absence of hot, young stars after a 
burst, and the latter is the case for an e(a)-type galaxy since the 
nebular emission lines are highly extinguished. Hence strong 
HS absorption equivalent width tends to signal the presence of 
starburst; high reddening further implies the starburst could 
be younger than what dust-insensitive indicators of star for- 
mation would suggest. 

Since the inferred timescales, < 1 Gyr, of starburst are 
shorter than cluster-crossing timescales (typically a few Gyr), 
the higher abundance of starbursting e(a) galaxies compared 
to the field fraction strongly suggests the connection between 
enhanced star-forming activity and infalUng cluster substruc- 
tures, as discussed in § 3. 

6.2. The Cause of Enhanced Star-Forming Activity 

Fig. 12 compares the cluster and field [O ll]-emitting sam- 
ples, showing the e(a) abundance in Abell 85 1 is roughly at 
the same level as z ~ 1 field galaxies among [O ii] emitters. 
Since the e(a) fraction significantly declines from z ~ 1 to 
z 0.1, it is likely that e(a) galaxies are fractionally more 
abundant in Abell 85 1 than in the field at a similar redshift. 
We argue that this enhancement of cluster e(a) population is 
Ukely caused at the expense of continuously star-forming e(c)- 
type galaxies. First, our constraint on the burst age, < 1 Gyr, 
of e(a) galaxies is consistent with the triggering of starburst 
upon cluster infall, whose timescales are on the order of a 
few Gyr. Although the timescale argument does not necessar- 
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ily prove the greater frequency of starbursts upon infall [e.g., 
upon infall, e(c) galaxies could drop out of the [O ll] -emitting 
sample at a faster rate via truncation of star formation than 
the e(a) production rate via starburst], a higher burst fraction 
is consistent with post-starburst (E+A) galaxies being more 
abundant within the total galaxy population in clusters than in 
the field at intermediate redshifts (Dressier et al. 1999; Tran 
et al. 2003, 2004, but see Balogh et al. (1999)). Second, star- 
burst galaxies almost certainly require gas-rich progenitors. 
The only class of galaxies with reasonably ample fuel supply 
is star-forming e(c) galaxies. Due to their short lifetime, e(b) 
galaxies would not make a significant contribution. Third, the 
e(c)-to-e(a) transformation appears consistent with their mor- 
phology. Poggianti et al. (1999) find e(a) and e(c) galaxies in 
z ~ 0.4 clusters tend to show late-type disk morphology, but 
e(a) galaxies appear much more likely to be associated with 
mergers or strong interactions. Balogh et al. (2005) also find 
field e(a) galaxies in z w 0. 1 to show disk-dominated morphol- 
ogy. Furthermore, their distributions in the [Oll]-H(5 plane 
(Fig. 1 1) indicate e(a) galaxies are neither rare or particularly 
special in their properties. As clearly seen in the SDSS field 
sample, at least some galaxies in the e(a) class simply con- 
stitute a high-H(5 extension of e(c) galaxies following an evo- 
lutionary track for a small- to moderate-£'(B-y), r ~ several 
Gyr exponentially-declining model, which is a typical SFH 
for spiral galaxies. Perhaps the e(a) galaxies of Balogh et 
al. (2005) constituted of a mixture of H(5-strong, but other- 
wise normal e(c)-type galaxies and real dusty disk galaxies 
observed during an episode of starburst. The set of evidence 
here naturally points to the underlying connections between 
e(a) and e(c) galaxies — it is reasonable that triggering of star- 
bursts in gas-rich e(c) galaxies gives rise to e(a) galaxies in 
these contexts. 

The process of e(c) to e(a) transformation must induce in- 
creased amount of dust in e(a) galaxies. At high redshift be- 
yond z > 0.7, mergers are presumably more frequent, and 
infrared-luminous sources appear to dominate star-forming 
activity (e.g., Le Floc'h et al. 2005), implying the importance 
of star-forming activity obscured by dust. These observations 
indicate the apparent connection between dynamically active 
environment and its effect on inducing dusty starburst. Al- 
though major mergers are not likely to result in disk systems, 
merger-induced dusty starburst has been proposed as a mech- 
anism for reproducing e(a) spectrum (Shioya & Bekki 2000); 
dust plays an essential role in prolonging the e(a) spectrum 
phase by extinguishing emission lines from star-forming re- 
gions. It is known that most vigorously star-forming galax- 
ies emit strongly in infrared through dust reprocessing of UV 
light from young stars. Thus it is not surprising that e(a) 
galaxies are dusty, if they are in fact starburst galaxies. Fur- 
thermore, the selectively higher extinction of light from very 
young stellar population (Poggianti et al. 2001) would tend to 
yield high reddening values, if the nebular lines originating 
from vigorous star-forming activity can be measured at all. 
Our spectral evolutionary tracks assumed constant reddening 
E(B-V) throughout the lifetime of a starburst galaxy. If, in- 
stead, the amount of extinction has a decreasing trend with 
the age of the stellar population sampled for observation, the 
fact that e(a) galaxies are dustier just reflects the significant 
presence of very young stellar populations on top of under- 
lying older stellar populations. Therefore the high reddening 
values favor a picture of e(a) galaxies as starbursts by infer- 
ring the presence of very young stellar populations but may 
not help us constrain one mechanism over others as to what 
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Fig. 18. — Spatial distribution of cluster [Oil] emitters in terms of their 
D99 spectral class: e(b) (blue squares), e(c) (green circles), e(a) (red dia- 
monds), k (brown inverted triangles), and k+a/a+k (magenta triangles). Black 
crosses denote [On] emitters for which a spectral class cannot be determined. 
Filled symbols denote the spectral class derived from better than 2a measure- 
ments of hne indices, whereas open symbols denote that from < 2(t measure- 
ments. Smaller, slightly faded crosses are MLFOO [O II] emission-Une can- 
didates not observed spectroscopically; these are g — i color-coded by blue 
(1 < g-i), green (1 < g — i < 2), red (g-i > 2), and gray (the uncertainty in 
g — i is > 0.5). The contours are the same as in Fig. 8. [See the electronic 
edition of the paper for a color version of this figure.} 

actually causes starbursts, i.e., dusty starbursts do not have to 
be caused by mergers or interactions. 

The question yet to be answered clearly is what physi- 
cal mechanisms are causing the cluster [Oil] emitters to go 
through starbursts. To address this, we see in Fig. 18 the spa- 
tial distribution of cluster [O ll] emitters with their D99 spec- 
tral class. A remarkable observation is the high degree of clus- 
tering of e(a) galaxies near or inside the infalling groups, i.e., 
the southern and north-western clumps of galaxy overdensi- 
ties, and the periphery of cluster central region. The appar- 
ent increase of galaxies showing starburst signatures in their 
spectra near the KOI threshold density (Fig. 3) offers another 
evidence for such a trend. In galaxy groups beyond the clus- 
ter virial region, both major mergers (Icke 1985; Mihos 1995; 
Bekki 1998) and galaxy harassment (Moore et al. 1999) are 
probable, but the effect of intracluster medium is expected to 
be less effective. Therefore galaxy-galaxy interactions offer a 
natural explanation for the starburst triggering mechanisms in 
those environments. The best way to confirm this hypothesis 
is to look at their morphology for signs of interactions. Some 
high-resolution Wide Field Planetary Camera 2 (WFPC2) im- 
ages from Hubble Space Telescope {HST) exist for Abell 851, 
but unfortunately the infalling substructures seen in Fig. 18 
are outside of their coverage. 

Turning our attention to the cluster central region, we see 
e(a) galaxies residing mostly in the periphery of the high- 
est galaxy-density contour (Fig. 18). Although some of these 
[O ll] emitters could be 2D-projected foreground/background 
galaxies in moderately dense infall group environments along 
the filament outside the core, the [Oil] emitters are mostly 
seen in the lowest velocity bin in Fig. 8, indicating their radial 
velocities are relatively close to the cluster systemic. In the 
core environment, both the intracluster gas and gravitational 
potential of the cluster can interact with cluster galaxies. Ma- 
jor galaxy-galaxy mergers are expected to be very rare, since 




Infalling Faint [O ii] Emitters in A851. II. 



17 



TABLE 4 

Cluster [O ii] Emitters with HST WFPC2 Morphology 



A iffT T7AA 


ARA 


AUbC 


D99 


S97 


Hubble 


D 


• 

Interpretation^ 


Comment 


UJ 


( ) 


( ) 


Class 


LU 


Type^ 








1448 


+0.76 


-3.94 


e(a) 


314 


Sbc 


2 


Tidal feature? 


Face-on disk or disturbed spheroid? Shells or arms? 


1640 


+1.32 


-3.53 


e(a) 


2022 


Sb 







Multi-armed disk galaxy; bright central disk, small bulge 


2095 


+1.55 


-2.48 


e(b) 


499 


Sc 


1 




Small, highly inclined disk; sUght structure visible 


2303 


+0.46 


-1.94 


e(b) 


N/A^ 










2978 


+O.70 


-0.61 


e(a) 


2016 


E 









3005 


+0.55 


-0.50 


e(a) 


620 


Sb/SO? 


1 


Tidal feature? 


Small, high SB spheroid; featuring disk w/arms, or tidal 


3215 


+0.06 


-0.21 


e(a) 


497 


Sbc 


2 


Merger 


Visually disturbed, high SB spiral; probably merger with tail at 1 1 


3568 


+1.31 


+0.63 


e(a) 


34 


SBab 







High SB bulge or nucleus; weak arms from bar 


3600 


+0.27 


+0.55 


k+a 


122 


Sd 


3 


Merger 


Disk system with tails; bright blobs — likely merger 


3604 


+0.35 


+0.67 


a+k 


36 


Sd 


2 


Merger 


Very high SB, disky object w/ 2 blobs and tail 



Note. — First four galaxies are in Field 2 of the SmaU et al. (1997) survey which sampled galaxies in the "outskirts" of AbeU 851, and the latter 
six are from the cluster central region from the same survey. 

" MLFOO object ID from their [O II] emission-line search. 
Relative to 9''42"'56!'2, 46°59' 12^0 (J2000). 

D99 spectral class derived from our measurements. The classification derived from < la measurements are italicized. 

Cross-identification with Small et al. (1997). The IDs from the pointing outside the cluster center, i.e., C10939+47 Field 2 in their catalog, are 
italicized. 

* The standard Hubble classification scheme (E, SO, Sa, Sb, etc.). 

* Disturbance class: 0, normal; 1, moderate asymmetry; 2, strong asymmetry; 3, moderate distortion; 4, strong distortion. 

8 Interpretation of disturbance classes, roughly divided into (M) merger, (1) tidal interaction with neighbor, (T) tidal feature, (C) chaotic, and (!) 
remarkable. 

No matching object found in the morphology catalog. 



the deep gravitational potential well of cluster means that the 
galaxies close to the cluster dynamical center are moving very 
fast, reducing the probability for a low-speed encounter with 
another galaxy. Yet galaxy harassment, in the sense of high- 
speed encounters between galaxies, would still be effective in 
the central region. The morphological analysis using HST has 
been done for many galaxies in the central region of Abell 85 1 
by well-trained morphologists. In Table 4 we list a common 
subset from the catalog by Smail et al. (1997) and our [Oil] 
emitters. Of particular interest here are those objects found 
within the cluster central region (the latter six in the list). Four 
out of six galaxies have their disturbance class greater than 
0, meaning there is some morphological sign of interaction 
in those galaxies. Of these, three are presumably mergers, 
although we are not sure whether the Small et al. interpre- 
tations distinguish between high-speed (i.e., harassment) and 
low-speed galaxy-galaxy encounters (i.e., major/minor merg- 
ers). As the comment indicates the remaining one is not a 
clear-cut case, and the cause of its disturbance is not clear. 

In any case, the existing morphological data in Abell 851, 
though limited in quantity, indicate a high fraction of the clus- 
ter e(a) galaxies are involved in interaction. As discussed ear- 
lier, Poggianti et al. (1999) already noted a higher incidence 
of mergers or interactions in e(a) galaxies compared to e(c) 
galaxies from the full Smail et al. morphology catalog of in- 
termediate redshift clusters. We further note that most of the 
[O II] emitters in the central region are morphologically clas- 
sified into disk galaxies, yet their spectral signature impUes a 
current or recent starburst. In the presence of high reddening, 
their optical colors and the strength of [Oil] emission may 
not give an obvious indication of starburst, so the apparently 
"normal" spiral galaxies in the cluster center, which are pho- 
tometrically more akin to field spirals than cluster spirals in 
nearby clusters (Andreon et al. 1997), could be going through 
bursts. 

Of course it is of no surprise that starbursts are induced by 
some form of interaction; galaxy groups provide ideal envi- 
rormients for merger-induced star-forming activity. We now 
speculate if the occurrence of starbursts could be enhanced 



specifically in cluster environments. X-ray observations show 
that Abell 851 is going through an on-going cluster-scale 
merger process, in which two X-ray emitting cores merg- 
ing in about a few hundred Myr (Schindler & Wambsganss 
1996; Schindler et al. 1998; De Filippis et al. 2003). Recent 
simulations have shown that galaxy clusters in ACDM cos- 
mology have two distinct phases in mass accretion history, 
merger-dominated and accretion-dominated phases (Zhao et 
al. 2003). Major cluster mergers, involving almost doubling 
their total mass, is not rare at the redshift of Abell 85 1 (Ta- 
sitsiomi et al. 2004). Inferred star-forming property of cluster 
galaxies would show differences in each phase, adding pos- 
sible complications to a pure infall scenario. The apparently 
complex spatial as well as velocity distributions of [O ll] emit- 
ters in Abell 85 1 perhaps reflects such cluster-scale mixing of 
galaxy populations. Applying dynamical methods to such a 
dynamicaUy active cluster poses difficulties as well; velocity 
dispersions are not a good indicator of cluster properties for 
dynamically active clusters going through mergers (Crone & 
Geller 1995).'° 

On the flip side of the complexity, there are some evidence 
that dynamicaUy active clusters host a larger number of star- 
forming galaxies and AGNs (Lotz et al. 2003; IVIiller & Owen 
2003; Rakos & Schombert 2005). As mentioned previously 
the existence of old galaxy population tracing clumpy dis- 
tribution of substructures along filaments is generally con- 
firmed in wide-field survey of Abell 851 using both photo- 
metric redshift method (KOI) and cluster red-sequence tech- 
nique (Kodama et al. 2005). The radial velocity distributions 
show weak, yet positive evidence for kinematic segregation 
between the cluster [On] emitters and D99 k-type galax- 
ies, when the samples are restricted within the core region 
(Fig. 3). Their spatial distributions, again restricted within 
the central region, clearly show an abundance of old galaxies 
where the cluster [0 11] emitters only sparsely exist. Hence it 

As a result, the estimates for the virial radius of Abell 85 1 vary widely, 
and its unique dynamical center can not be rigorously determined. This is a 
major reason why we have stayed clear of analysis using cluster-centric radii. 
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is natural to interpret many [O ll] emitters to be relatively new, 
recently accreted galaxies onto the cluster potential which al- 
ready hosts many evolved galaxies. The non-Gaussian ve- 
locity distribution of cluster [O ll] emitters also provides evi- 
dence that they represent a population that has not been virial- 
ized in the cluster potential. In general, it is unlikely that red, 
gas-poor galaxies such as cluster ellipticals get rekindled in 
starbursts, but a dominant contribution to cluster star-forming 
activity can easily come from gas-rich spiral galaxies infalling 
onto a cluster Hence, the enhanced star-forming activity in 
the infall region, like that seen in the southern substructures in 
Fig. 18, probably is due to an increased frequency of galaxy- 
galaxy interactions between quiescent yet gas-rich galaxies, 
as the infalling groups themselves interact with nearby groups 
on the way to the cluster center This is consistent with the 
high abundance of [Oil] emitters between the apparent in- 
falling groups in the southern filament. 

Yet clusters offer means of interactions other than the ones 
between galaxies — the dense cluster gas in the core could 
ram-strip galaxies of their fuel supply, and the steep clus- 
ter gravitational potential could tidally interact with infalling 
galaxies. In CL0024H-1652 at z w 0.4, Moran et al. (2005) 
found that morphologically early-type galaxies with large 
mass showed enhanced Balmer absorption lines within the 
cluster virial radius, while some fainter counterparts close 
to the virial radius actually showed [O ll] emission. These 
authors rule out merger-induced star formation after inspec- 
tion of their morphology (though they mention harassment 
is a potential triggering mechanism). In local Coma cluster, 
Poggianti et al. (2004) recently found an abundance of young 
post-starburst galaxies near the edges of X-ray substructures, 
indicating the interaction between cluster galaxies and intra- 
cluster medium could trigger star-forming activity. These ob- 
servations imply the importance of interactions between in- 
falling galaxies and environments specific to clusters. The 
X-ray observation by De Filippis et al. (2003) was probably 
not deep enough to trace the cluster gas much beyond the very 
center of Abell 851 (Fig. 19). 

To summarize, in view of highly dynamically active state 
of Abell 851, the enhancement of star-forming activity seen 
in this cluster may be the result of accelerated gas consump- 
tion within infalling gas-rich galaxies driven by the galaxy- 
galaxy interactions in infalling group environments. As a fa- 
vored triggering mechanism of starbursts, galaxy-galaxy ha- 
rassment, as opposed to a major merger, appears consistent 
with preserving the disk morphology of cluster [O ll] emit- 
ters, although the direct confirmation with our [O ll] -emitting 
sample is limited only to the cluster central region due to the 
limited coverage of HST observations. For the rest with sub- 
tle or no morphological indication of interaction in the clus- 
ter central region, we can only speculate if an increased in- 
homogeneity in the distribution of intracluster gas/potential 
makes triggering star-forming activity more effective, such as 
in shocking the interstellar medium of galaxies. 

6.3. Possible Role ofAGNs as Firemen 

Some optically-luminous [O ll] emitters in the cluster cen- 
tral region shows a glimpse of the fate of cluster star-forming 
galaxies. To see how the [O ll] emitters are distributed relative 
to the intracluster gas, in Fig. 19 we overlay X-ray emission 
contours, kindly provided by B. De Filippis. We see very few 
[O II] emitters in the region where the X-ray emission is very 
strong, although galaxies of the types associated with old stel- 
lar population are abundant (see Fig. 8). More remarkably. 
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Fig. 19. — The disti'ibution of cluster [Oil] emitters in the central region 
of Abell 851 with the HST WFPC2 images. The north is up, and the east is 
in the left. The yellow cross denotes 9'"42""56?2, 46°59'12"0 (J2000). from 
which ARA and ADEC are measured in Fig. 8. The overlayed contours in 
cyan are X-ray emission map from De Filippis et al. (2003). Those galaxies 
identified to be AGN-like are marked by black squares; otherwise the colors 
of the symbols denote D99 spectral classes as in Fig. 18. [See the electronic 
edition of the paper for a color version of this figure.} 

the [O II] emitters classified as AGN-like (in black squares) 
appear to align themselves in a preferred direction from the 
south-west to the north-east. Whether this preferred axis has 
any physical significance is not obvious both from the figure 
and from the X-ray temperature map presented by De Filip- 
pis et al. (2003). Nonetheless we note that the cluster central 
field presented in Fig. 19 is a small fraction of the total survey 
area for our [O ll] emission search but hosts six of the seven 
AGN-like galaxies in our sample; in comparison, only 20 out 
of 111 confirmed cluster [Oil] emitters and 186 out of 808 
KOI cluster members are in the region shown in Fig. 19. 

We reiterate that our AGN identification is not very com- 
plete, i.e., the line-ratio diagnostics could be done only for 
a subset of all the cluster [O ll] emitters. Furthermore, none 
of the AGN-like galaxies are spatially coincident with X-ray 
point sources identified by De Filippis et al. (2003); only one 
g-i> 2{0 II] emission candidate (which was not followed up 
spectroscopically) from the MLFOO narrowband survey coin- 
cided with the point source #8 in De Filippis et al. The op- 
tical diagnostics are often inadequate for secure AGN iden- 
tification; Martini et al. (2002), for example, have reported 
in Abell 2104 the optically-selected AGNs fraction is by a 
factor of 5 lower than that of the X-ray-selected — some X- 
ray-selected AGNs do not even show strong emission lines. 
On the other hand, a high X-ray column can give rise to X- 
ray-obscured, optically-identified AGNs (Brandt et al. 2004). 
Based on the fluxes of the X-ray point sources detected in 
their imaging, the XMM exposure time of 5 1 ks should have 
easily detected AGNs at the flux level of /z(0.3 - 10.0 keV) > 
10" erg s~' cm"^. To our knowledge it is not clear whether 
a reasonable level of high X-ray column can obscure AGNs 
to put our galaxies below this level of detection. Nonetheless 
the optical line ratio diagnostics does suggest a possibility that 
these AGN-like galaxies could harbor low-luminosity AGNs. 

If these are in fact AGNs, a possible clue may be provided 
by a recent study by Yan et al. (2005), reporting that the [O ll] 
emitters in red sequence and post-starburst galaxies in their 
low-z SDSS sample seem to host AGNs/LINERs. In an in- 
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dependent study, Tremonti et al. (2006, private communica- 
tion) similarly finds signs of AGN activity in their SDSS post- 
starburst galaxy sample. Spectral decomposition of an SDSS 
AGN sample by Vanden Berk et al. (2006) showed post- 
starburst activities to be much more common in broad-line 
AGNs. Interestingly, three of the AGN-like galaxies in our 
sample are classified into k-na/a-nk, a post-starburst type, in the 
D99 spectral classification scheme, although the fact that they 
made into our [O ii] -selected sample means they have residual 
[O II] emission. We also find most of the AGN-like galaxies 
in our sample reside in local environments whose galaxy den- 
sities are higher than the KOI threshold above which cluster 
galaxy colors are typical of red, non-star-forming population 
(although we must point out the KOI threshold characterizes 
the dense cluster central region, as well as infalling substruc- 
tures where we find little evidence for AGN-like galaxies). 
Furthermore, their location on the color-magnitude diagram 
indicates they are among the most luminous [O II] emitters, 
but their g - i color is either indistinguishable from the rest of 
the [O ll]-emitting population or consistent with being on the 
red sequence (Fig. 6). With the indications that the cluster red 
sequence building up from the luminous end in a downsiz- 
ing fashion (De Lucia et al. 2004; Tanaka et al. 2005), there 
is a curious connection between AGNs as a mechanism for 
shutting off star formation and the property of host galax- 
ies. Hence the distribution of AGN-like [O ii] emitters in the 
color-magnitude diagram, along with their relatively smaller 
[O II] emission equivalent widths (Paper 1), is consistent with 
the scenario in which the star-forming activities within lumi- 
nous galaxies are shut off by AGNs. Similar observations of 
cluster star-forming galaxies over a wide range of redshifts 
could provide evidence for or against the role of AGNs in 
causing the downsizing trend, where quenching of star for- 
mation progressively move toward lower luminosity systems 
at lower redshifts. 

Yet, it is a mystery as to why AGN-like galaxies are not 
seen in the [O ll] emitters outside the cluster central region. 
The KOI threshold effect, if real, implies that galaxies can 
transform from star-forming to non-star-forming in group en- 
vironments in infall regions. Six out of seven AGN-Uke galax- 
ies in our sample are in high-density environments above the 
KOI threshold but are preferentially located in the cluster cen- 
tral region. Ruderman & Ebeling (2005) have generally found 
an excess of X-ray-selected AGNs in the core region in both 
relaxed and morphologically-disturbed clusters (plus another 
milder peak near virial radius in relaxed clusters). These ob- 
servations might suggest the significance of the locations of 
galaxies within clusters, as well as the effect of local environ- 
ments, in transforming the star-forming properties of cluster 
galaxies. 

There also is an ambiguity as to whether their broadband 
/ magnitudes are truly representative of the host galaxies. If 
AGNs are active, it is likely that they make a significant con- 
tribution to the continuum light in the optical, making host 
galaxies brighter than they really are. That the most AGN-Uke 
galaxies in our sample are among the brightest [O ii] emitters 
might simply be a result of additional contributions to the con- 
tinuum from AGN activities. 

In view of the general difficulty in systematic AGN selec- 
tion, the AGN-like galaxies in our [O ii]-emitting sample do 
not offer direct evidence for the connection between quench- 
ing of star formation and AGN activity. Yet, the observed 
trends discussed here are all consistent with AGNs playing 
some role in the production of bright cluster ellipticals with 



little star-forming activity. In any case, a more robust AGN 
identification and follow-up studies of these objects would be 
necessary to confirm this hypothesis. 

6.4. Interpretation of Faint [Oil] Emitters 

The analysis of composite spectra based upon sampling on 
the ;-S plane has shown relatively little variation with local 
galaxy density, except for the faintest [O II] emitters residing 
in the low-density environment, i.e., the LF subsample. The 
[O II] emitters in the other subsamples are mostly consistent 
with the e(a)-type and typically show higher H(5 absorption 
equivalent widths at ^ 4-5 A, with relatively little scatter. 
The faint [Oil] emitters in low-density environments show 
the [O II] and H5 equivalent widths typical of e(b) spectrum 
indicating a very young, pure starburst. 

A^-body simulations by Moore et al. (1999) demonstrated 
that the response of a spiral galaxy to the tidal disruption 
from the cluster potential and galaxy-galaxy interactions — 
galaxy harassment — varies depending on the concentration of 
the mass distribution and the disk scale length. Low surface 
brightness galaxies respond dramatically to the influence of 
both the local and global variations in the gravitational po- 
tential, while high surface brightness galaxies remain stable. 
We lack morphological information to see how ; magnitude 
and mass concentration are related systematically in our sam- 
ple, though faint g-i < 1 [O ll] emitters generally have small 
isophotal area, implying their dwarf identity (MLFOO). The 
relative depletion of HF [O II] emitters in Fig. 5 could be evi- 
dence for destruction of star-forming dwarf galaxies in dense 
regions due to galaxy harassment. On the other hand, such an 
absence of faint star-forming population does not have to indi- 
cate the total absence of galaxies in the particular subset. The 
depletion of faint [O II] emitters in dense environments may 
just be a reflection of highly suppressed star-forming activity 
in dwarf galaxies. The dense intracluster medium can pro- 
vide effective gas-removal mechanisms, such as ram-pressure 
stripping, which is also a very likely scenario for dwarfs with 
a more diffuse matter distribution. In this case galaxies are 
not destroyed but could simply fade out of i band flux limit as 
star formation ceases. 

The contribution of dwarf galaxies to the total galaxy popu- 
lation has been studied in terms of galaxy luminosity func- 
tions (e.g., de Propris et al. 1998a; de Propris & Pritchet 
1998b; Trentham & Tufly 2002; Valotto et al. 2004; Popesso 
et al. 2005, 2006). Up to and below the magnitude Umit rel- 
evant for the [O ii] emitters studied in this paper, there are 
evidence for a steeper faint-end slope in nearby clusters than 
that of field luminosity functions, and a higher abundance of 
dwarf galaxies in clusters has been proposed as the cause of 
such steep cluster luminosity functions (e.g., de Propris et al. 
1998a; de Propris & Pritchet 1998b). The detection of such 
a trend seems not universal, however, and there is even a 
suggestion that sample selection biases could cause a steeper 
faint-end slope (Valotto et al. 2004). Although it is still not 
clear how the general form of faint-end luminosity function 
differs systematically between clusters and field regions, type- 
specific luminosity functions appear to depend strongly on en- 
vironment regardless of the precise form of luminosity func- 
tion for the entire population: The brighter end is dominated 
by non-star-forming population, while the faint-end is very 
sensitive to the relative abundance of both star-forming and 
non-star-forming populations (Goto et al. 2002; De Propris et 
al. 2003). 

Two major surveys, SDSS (Goto et al. 2002) and 2dFGRS 
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(De Propris et al. 2003), disagree in their measurement of the 
faint-end slope of the cluster luminosity function — flatter in 
SDSS and only sUghtly steeper in 2dFGRS. As far as the in- 
terpretation of where the star-forming dwarfs in high-density 
environments go, the behavior of the faint-end slope makes 
a profound difference. If the disappearance of star-forming 
faint dwarfs is simply due to their fading as in stripping of gas 
supply, an increase in the number of non-star-forming dwarfs 
would match the decrease in the star-forming population, and 
the shape of overall luminosity function would not change 
as much, consistent with 2dFGRS, finding non-star-forming 
population to increase and make the faint-end slope steeper 
(less complications arising from flux limits being sensitive to 
galaxies with different stellar populations). On the other hand, 
galaxy harassment can either destroy galaxies as in tidal de- 
bris (Moore et al. 1999) or create even smaller galaxies as in 
tidal dwarf galaxies (Popesso et al. 2006). Hence the behavior 
of the faint-end slope of luminosity function alone would still 
not distinguish these dwarf galaxy production models. Our 
sample is also insensitive to the surviving dwarf galaxies in 
the cluster central region, since they are likely to be stripped 
of their gas supply (Mori & Burkert 2000) and show no star 
formation. 

Yet our data could offer insight as to at which point of 
cluster formation history the abundant cluster dwarf popula- 
tion is really produced. Since the metal enrichment ceases in 
such gas-poor dwarfs, comparing the gas-phase metallicities 
of dwarf galaxies in Abell 85 1 and local clusters may tell us 
if gas-rich dwarfs in z = 0.4 clusters could be the progenitors 
of the abundant gas-poor dwarf galaxies in local clusters. We 
wish to present the metallicity analysis in our forth-coming 
paper of the series. 

In short, our finding in Abell 851 merely shows high- 
density environments to be hostile to dwarf starburst galax- 
ies. The reported increase in non-star-forming dwarfs in clus- 
ter central region (Thompson & Gregory 1993; De Propris 
et al. 2003; Popesso et al. 2006) does not necessarily prefer 
one interpretation over another as to the fate of infalling star- 
forming dwarfs. The forth-coming metallicity analysis may 
provide some clues for addressing the progenitor-offspring is- 
sues surrounding cluster dwarf galaxies. 

7. SUMMARY 

We have reported on a spectroscopic follow-up study of 
narrowband-selected [On] emitters in Abell 851 from the 
MLFOO catalog, focusing on their local environments, kine- 
matics, and SFHs. Along with similar studies (e.g., Balogh 
et al. 2004a,b; Kodama et al. 2004; Rines et al. 2005), the 
cluster [On] emitters reaffirm the fundamental importance 
of local envirormients: The fraction of galaxies with [O ii] 
emission is a strong function of local galaxy density where 
few [O II] emitters are found in dense environments, yet the 
strength of [O II] equivalent width does not vary strongly with 
local galaxy density. The trends are consistent with the mech- 
anisms that halt star-forming activity on short timescales or 
are sufficiently rare. In high-density environments, the cluster 
[O II] emitters appear to show distinct behaviors depending on 
their luminosity: For optically-luminous [O II] emitters in the 
cluster central region, there is some weak evidence for AGNs 
playing a role in suppressing star formation in the host galax- 
ies. On the fainter end, we note the lack of faint [O ii] emit- 
ters in dense environments; the destruction of dwarf galaxies 
in galaxy harassment (Moore et al. 1999) and/or simple fad- 
ing of emission Unes due to quenching of star formation in 



ram-stripped dwarfs (Mori & Burkert 2000) can cause such a 
deficit. The luminosity-dependent behaviors of star-forming 
galaxies strongly suggest normal disk galaxies, which are 
predominantly e(a) galaxies, and dwarf galaxies, which are 
mostly pure bursting e(b) galaxies, respond quite differently 
within dense cluster environments. 

The cluster [O ll] emitters have shown the kinematic prop- 
erties consistent with a non-virialized population; their veloc- 
ity distribution deviates significantly from a Gaussian. The 
skewed distribution toward higher velocities is also consistent 
with a high fraction of the [O II] emitters being at first infall, 
as opposed to a backsplash population which may have gone 
through the cluster central region in the past (e.g.. Gill et al. 

2005) . We observe a greater fraction of cluster [Oil] emit- 
ters with high radial velocities, compared to the emission-Une 
galaxies in local clusters (Rines et al. 2005; Pimbblet et al. 

2006) . In addition, their spatial distribution is highly asym- 
metrical with respect to the cluster center. These provide evi- 
dence for a signification fraction of the cluster [O ll] emitters 
constituting an infalUng population. 

The spectral evolution modeling and the distribution of the 
cluster [O ll] emitters in the [O Il]-H(5 plane highlight a dom- 
inant presence of dusty starburst galaxies [i.e., e(a)-type in 
the D99 scheme]. A high abundance of e(a) galaxies in clus- 
ters was already remarked by Poggianti et al. (1999), whose 
subsequent studies identified them as dusty starbursts, pri- 
marily caused by selective extinction of the light from very 
young stellar population (e.g., Poggianti & Wu 2000; Pog- 
gianti et al. 2001). A combination of stellar population syn- 
thesis and nebular photoionization codes favors models with 
star-forming timescales of order < 1 Gyr for dusty galaxies, 
E(B-V) > 0.5, for covering the distribution of cluster [O ll] 
emitters in the [O Il]-H(5 equivalent width plane. The range of 
timescales, along with a higher occurrence of e(a) galaxies in 
Abell 85 1 compared to the field, is consistent with scenarios 
in which triggering of the starbursts is cluster-related. Exist- 
ing morphological studies (Small et al. 1997; Poggianti et al. 
1999) further provide evidence for e(a) galaxies arising at the 
expense of regular star-forming galaxies [i.e., e(c)-type] by 
preserving their disk-hke morphology. The redshift trend of 
the relative fractions of [O ii] emitters in various D99 spec- 
tral classes also suggest a higher incidence of e(a) galaxies in 
Abell 85 1 compared to the field at a similar redshift. 

The enhanced star-forming activities observed in Abell 85 1 
may have some connections to the highly dynamically active 
state of the cluster. We specifically note a high abundance 
of [O II] emitters between two infalling groups located in the 
southern filament of Abell 85 1 . The proximity of two groups 
during the process of an accretion onto the cluster may have 
caused the galaxies to interact at a higher frequency, although 
we lack high-resolution morphology to confirm signatures of 
interaction in these galaxies. An inhomogeneous intraclus- 
ter gas distribution (De Filippis et al. 2003) might also in- 
crease the frequency of triggering of star-forming activity in 
the galaxies in the cluster central region; a direct confirmation 
of this additional hypothesis would need a larger inventory 
of morphology and spectra for the galaxies spatially coinci- 
dent with a deep X-ray emission map. On the other hand, the 
archival HST morphology suggests galaxy-galaxy harassment 
is very effective in inducing a starburst even in the galaxies 
residing within the cluster central region. Galaxy harassment 
is also preferred over major mergers because of the assumed 
preservation of disk structures. 
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APPENDIX 

SPECTRAL EVOLUTION MODELING 

We synthesized model galaxy spectra incorporating differential extinction of stellar continuum and nebular emission compo- 
nents. We used a combination of a stellar population synthesis code, GALAXEV version 2003 (Bruzual & Chariot 2003), and a 
nebular photoionization code, version 96.01a of Cloudy (Ferland et al. 1998). First, for each predefined SFH, GALAXEV was 
used to synthesize a series of spectra at various ages. Three SFHs are explored: instantaneous burst; constant burst of finite 
duration t, after which there is no star formation (i.e., truncation models); and exponentially-declining SFHs with various r (no 
gas recycling). The analytical forms of these SFHs are described in Bruzual & Chariot (2003). For the latter two SFHs, we used 
the values of r = 10, 10^, 10^, and 2 x 10^ Myr. We used the GALAXEV models with the Padova 1994 evolutionary tracks 
and assumed a Salpeter initial mass function {M = 0.1-IOOMq) at solar abundance {X = 0.70, Y = 0.28, and Z = 0.02). The 
standard GALAXEV output gives the production rates of hydrogen-ionizing photons, 2(H), which are scaled to the total mass 
of an ionizing source (log7V4*/A^© ~ 2-9). The mass of ionizing source therefore enters only as a scaling factor for Q(H). 
Thus constructed ionizing source spectrum energy distribution and properly scaled (2(H) were fed to Cloudy for photoionization 
simulations. 

In Cloudy, our model H II region is a single hydrogen-ionizing source surrounded by a spherical distribution of constant-density, 
ionization bounded nebula at the solar abundance. We used the solar abundance and depletion factors assumed by Kewley et al. 
(2001). The hydrogen density (wh = 10 cm~^), the volume filling factor (e = 1), the covering factor (C = 1), and the inner radius 
(rin = 10"^ pc) to the nebula were all kept constant throughout the simulations with Cloudy. These parameters can be chosen 
rather arbitrary, since the ionization state of nebula is uniquely determined by the combination of the spectral energy distribution 
of ionizing source and the ionization parameter defined as 

6(H) 

where is the Stromgren radius and c is the speed of light, so long as the density remains below the Umit for coUisionally 
de-excitation coohng (see, e.g., Stasinska & Leitherer 1996). 

Dust modeling features were generally turned off within these codes. We manually attenuated the stellar and nebular compo- 
nents by a different amount of extinction, adopting the attenuation curve of Calzetti et al. (2000, Eq. (8)). It has been shown 
empirically that the stellar continuum suffers roughly half of the extinction suffered by the nebular gas, EiB-V)^ ~ 0.44EiB-V)g 
(see Calzetti et al. 1994, and reference therein). This ratio was varied [E(B-V)^JE(B-V)g = E^/E^ = 0.22, 0.44, and 0.881 to 
see the effect of differential extinction on our measurements of equivalent widths. This is motivated by the interpretation of e(a) 
spectrum by Poggianti et al. (2001), requiring selective extinction of light from stellar populations of various ages. A range of 
extinction E(B-V) = 0.0-1 .5 was explored. Note that E(B-V) = E(B-V)„ is the amount of extinction in the nebular component, 
being derived from a Balmer decrement (Paper I). Resulting synthesized spectra are used to measure [O ii] and H(5 equivalent 
widths using the flux-summing method (Paper 1). The advantage of our method presented here is that both emission and absorp- 
tion lines, each suffering a different amount of extinction, are modeled into the final spectra, which are measured with the same 
tool used for measuring other spectra. 

We checked our method of simulation by running a simple model that closely resembles the solar abundance, Mj = 2 x IQ'^Mq 
model published by Garcia- Vargas et al. (1995, Table 16), which has also been used by Moy et al. (2001) to check their simulation 
for consistency. Table Al shows that the results generally agree well with each other. 

Although our interest mainly lies in the evolution of spectral indices in the [O Il]-H5 plane, we need some constraint on the 
relevant range of ionization parameter u so that our models roughly resemble the ionization state of nebulae in the observed 
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TABLE A 1 

Comparison of Coupled-Model Simulations 



Models 



Line 


This Paper* 


MOl'' 


G95'= 


logH/3[ergs-'] 


38.79 


38.83 


38.84 


[Oll]A3727/H,9 


3.12 


3.10 


3.07 


[Oll]A5007/H/3 


0.32 


0.31 


0.25 


[Ol]A63()0/H/3 


0.07 


0.04 


0.04 


[NlI]A6584/H/3 


1.29 


1.34 


1.34 


[Sll]A6716/H/3 


0.70 


0.57 


0.58 


[Sll]A6731/H/3 


0.49 


0.39 


0.40 


[Slll]A9()69/H/3 


0.31 


0.36 


0.34 



Note. — The comparison of Z = Zq instantaneous 
burst model at 1 Myr and log m ~ -3 . 11 . The solar abim- 
dance is as defined in Garcia- Vargas et al. (1995). 



^GALAXEV (version 2003) & Cloudy (version 
96.01a). 

''PEGASE.2 & Cloudy (version 90.04) by Moy et al. 
(2001). 

■^Cloudy (version unknown) and see Garcia- Vargas et 
al. (1995) for their stellar population synthesis models. 
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Fig. Al . — Emission Une flux ratio [N ll]A6583/Ha vs. [O ll]A3727/[0 lll]A5007 for three evolutionary model tracks. The numbers along each track indicate 
the ionization parameter u for each model. For each evolutionary track, data points appear at 1, 10, 10^, 10^, and 5 x 10' Myr from high to low u. Gray data 
points indicate the cluster [O II] emitters whose line ratios are computed from their highest-S/N spectrum (filled dots) and secondary spectrum (crosses). Squares 

photons 2(H), and electron temperature T^, namely 



1/3 



where A(7;) « 2.8 x lO^^'^ClOV?;)^/^ (Stasinska & Leitherer 1996) and the Te dependence comes from Rs (which depends on the 
recombination coefficient). Since the temperature varies within a nebula, we simply take the value of u computed by Cloudy at 
the end of each simulation. Hence adjusting u is done through one free parameter, the mass of central ionizing source with which 
(2(H) scales proportionally. Since the strength of [Oii]A3727 emission depends on the relative abundance of different oxygen 
ionic species, we constrained this by comparing the ratio [O lll]A5007/[O ll]A3727 between the [Oil] emitting sample and our 
simulated models. In Fig. Al we see how solar-abundance models compare to observed galaxies. We first notice that any single 
model track between the ages f = (1-5) x 10^ Myr does not span the range of observed [O ll]/[Olll] entirely, suggesting a need for 
variety of SFHs to explain the wide range of observed hne ratios. Moreover, our models yield generally higher [N ii] A6583 /Ha. 
The ratio [Nil] /Ha is sensitive both to the relative abundance of N to all other species and to the overall metal abundance of 
nebula, but it has relatively little effect on the observables in the [Oll]-H(5 plane. Although the line ratio diagram implies the 
cluster [O ii] emitters may be at subsolar abundance in general, we chose solar-metallicity, M. = 4Mq models at E{B-V) = 0.5 
and E^/Eg = 0.44 as our fiducial set of models. 
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